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IMPLEMENTATION AND FLIGHT TESTS FOR DIALS 


INTRODUCTION 

In the early 70' s, LaRC's Terminal Configured Vehicle (TCV) program was 
formulated with an overall goal of improving the efficiency and capabilities of 
transport operations (ref. 1). A specific goal, in support of its broader objective, 
was to increase the capabilities of autoland systems. To achieve this goal, a 
research effort was begun to design, develop, and flight test an advanced autoland 
system which eventually became known as the Digital Integrated Landing System 
(DIALS). This is part one of a report which describes the implementation and flight 
tests and presents flight test data of the DIALS. 

Prior to the design phase, guidelines were formulated to establish the general 
system architecture for the DIALS and to take advantage of emerging technologies. 

The first guideline was to design the system using direct-digital design methods. 

Such a design could be implemented in the digital flight computers (which were to be 
used in future aircraft) without having to transform the control laws into the 
digital or discrete domain as is necessary for analog designs. A second guideline 
was to design an integrated controls system using modern control theory methods as 
opposed to the classical design methods currently used by the industry. Modern con- 
trol theory methods provide a means for designing multiple-input, multi-output 
systems whose control commands are coordinated to assist each other during a flight 
path maneuver. On the other hand, classical methods are restricted to single-input, 
single-output methods. A typical problem of classical designs is that the individual 
controls are not necessarily coordinated and can work against each other to reduce 
the overall closed-loop system stability and performance and to introduce unnecessary 
control activity. Another reason for using modern control methods was that these 
methods had evolved (over a period of 20 years or more) to a level where practical 
designs were considered to be achievable. However, a successful flight test of such 
a system had not been achieved. Thus the design, development, and successful flight 
test of a modern control design was necessary to establish this method as a useful 
tool for integrated control designs. The last guideline was to use the digital navi- 
gation signals of the Microwave Landing System (MLS) in the design rather than the 
analog signals of the current Instrument Landing System (ILS). The MLS, being 
developed by the Federal Aviation Administration ( FAA) , will provide signals with 
more accuracy and more volume coverage (ref. 2) than that provided by the ILS. Thus, 
the MLS was to be used to provide better performance and to allow for more flexi- 
bility in the final approach-to-landing path. 

The design and development of DIALS is covered in references 3, 4, 5, and 6. 

The basic equations with a general discussion of the DIALS control modes and flight 
tests with some flight test data are presented in reference 7. The details of the 
implementation of the DIALS in the flight control computers are not included. Thus, 
part of this report is to provide the implementation details of the DIALS in the 
flight control computers and present a more comprehensive selection of the fliqht 
test data. Part two of the report presents the data from the last five flight tests 
of the DIALS (those free of major software errors) and specifies the differences in 
the software system from flight to flight and run to run. 

In this part of the report a general description of the DIALS and the test 
vehicle will be given. Also, the sensor measurements that were used by DIALS will be 
described. A discussion of the equations implemented and the modifications of these, 
which were made during the flight hardware/software implementation process, follows. 


Next, a description of the test 
the results of the flight tests 
concluding remarks. 


procedures for a test run will he given, 
will be presented and discussed followed 


Finally, 
by some 


SYMBOLS LIST 
Measurements 


1 . Sensors for Control System 
$ roll attitude, rad 

0 pitch attitude, rad 


P 

q 

r 

EPR1 

EPR2 

h B 


B 

Th 

Th. 


RF 


LF 


es 


trim 


yaw or true heading, rad 

body-mounted acceleration specific force along x-body axis, ft/sec 

body-mounted acceleration specific force along y-body axis, ft/sec 

body-mounted acceleration specific force along z-body axis, ft/sec 

roll rate about aircraft x-body axis, rad/sec 

pitch rate about aircraft y-body axis, rad/sec 

yaw rate about aircraft z-body axis, rad/sec 

engine pressure ratio of the number 1 engine 

engine pressure ratio of the number 2 engine 

barometric altitude, ft 

radar altitude, ft 

barometric altitude rate, ft/sec 

position of the right forward throttle, deg 

position of the left forward throttle, deg 

stabilizer position, rad 

elevator servo position, rad 

right aileron servo position, rad 

left aileron servo position, rad 

aileron trim position, rad 


i 


MLS azimuth angle, deg 


MLS elevation angle, deg 


MLS range, ft 


calibrated airspeed, ft/sec 
rudder servo position, rad 


R. . 
trim 


rudder trim position, rad 


2. Processed 


acceleration along the x-body axis, ft/sec^ 

, o 

acceleration along the z-body axis, ft/sec 

2 

acceleration along the x-stability axis, ft/sec 

2 

acceleration along the z-stability axis, ft/sec 
total thrust of both engines, Klbs 

adjustment term to account for mean-sea-level attitude of runway and 
the Earth’s curvature, ft 

x-coordinate of the aircraft c.g. in the MLS coordinate frame, ft 

y-coordinate of the aircraft c*g. in the MLS coordinate frame, ft 

z-coordinate of the aircraft c.g. in the MLS coordinate frame, ft 

vertical velocity estimate from the second-order complementary filter 
or test vehicle, ft/sec 

vertical velocity estimate from the third order complementary filter of 
test vehicle (see appendix B) , ft/sec 

x-ciirdubate if MLS ground elevation antenna in MLS coordinate frame, ft 

y-ciirdubate if MLS ground elevation antenna in MLS coordinate frame, ft 

z-ciirdubate if MLS ground elevation antenna in MLS coordinate frame, ft 

x-coordinate of MLS ground DME antenna in the MLS coordinate frame, ft 

y-coordinate of MLS ground DME antenna in the MLS coordinate frame, ft 


z-MLS coordinate computed from radar altitude measurement, ft 
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3. Constants 


h R1 constant to incorporate the radar altimeter measurement into MLS processing 

non-dimenional H ' sS,in g, 

*GPIP x-coordinate of glidepath intercept point in the runway coordinate frame, ft 

H TDC distance from aircraft center-of-gravity to bottom of main landing gear 

wheels with gear struts extended, ft 

2KN2 constant to account for Earth’s curvature, 1/f t 

KZO constant for mean-sea-level atitude of runway, ft 

UOI constant equal to 1 divided by desired calibrated airspeed, sec/ft 

R o MLS range measurement processing constant, ft 

R c length of MLS c-band antenna cable to MLS receiver, ft 

^R/W true heading of runway, rad 

— — Additiona l Sen s ors and Variables for Plotting a nd Data Analysis 

a n specific force of normal accelerometer equal to the negative of f z , ft/sec 2 

a x measured specific force of longitudinal accelerometer equal to f ft/sec 2 

EPR aV lZT, ° f th t engine P^ssure ratio measurements from each engine (EPR1 and 

EPR2;, non-dimensional 

^ measured aerodynamic sideslip, deg 

6 e elevator surface position relative to stabilizer chord plane, deg 

6th th< fTh Vera9 a °l tl ? e left and right forward flight deck throttle positions 

' ^rf snd Th-^p ) , deg 

6 Th EST estimate of the throttle position equal to z LQ + Th Q , deg 

6 Th AFT the measured position of the aft flight deck throttle position, deg 

Ah the estimated altitude deviation from the desired glideslope equal to 

the negative of e Lx6 U 0 , ft 

Ay estimated distance from runway centerline equal to x f U ft 

6 o 9 L u 

Aip yaw deviation from runway centerline equal to iji - tjj deq 

^EST estimated yaw deviation from runway centerline equal to x 2 (57.3), deg 

A<P tk aircraft ground track angle deviation from runway centerline egual to 

measured track angle from inertial navigation system minus T , deg 

R / W r ^ 
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Filter Variables and Parameters 


1 . Longitudinal 


b 


L31 



constant gain on elevator servo position for vertical acceleration 
innovation, 1 /sec 

the i^ row j ^ column filter gain element for the aircraft state estimation 


X L1 

X L2 


X L3 

X L4 

X L5 

X L6 

*L7 

X L8 

X L9 

f 

LW. . 
ID 

W L1 


W L2 

W L3 

W L4 

W L5 


w 


L6 


w. 


L7 


LI 


L3 


L4 


L5 


L6 


L7 


perturbed pitch estimate, rad 

perturbed normalized inertial speed component along the the x-stability 
axis, non-dimensional 

perturbed inertial alpha, rad 

perturbed pitch rate, rad/sec 

normalized inertial runway x-axis component, sec 
normalized inertial runway z-axis component, sec 
perturbed thrust, Klbs 
perturbed throttle position, deg 
perturbed stabilizer position, rad 

the i th row j ^ column filter gain element for wind state estimation 
normalized wind gust component along the z-body axis, non-dimensional 
wind gust state, non-dimensional 
wind gust state, non-dimensional 

normalized wind gust compinent along the x-body axis, non-dimensional 

normalized steady state wind along the x-body axis, non-dimensional 

normalized steady state wind along the z-body axis, non-dimensional 

normalized shear wind component along the x-body axis, sec" 1 

pitch innovation, rad 

x position innovation, sec 

z position innovation, sec 

barometric altitude innovation, sec 

vertical velocity innovation, non-dimensional 

vertical acceleration innovation, sec"* 1 
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V L8 airspeed innovation, non-dimensional 

V L9 longitudinal acceleration innovation, sec -1 

y L l incremental pitch measurement, rad 

7^2 pitch rate measurement, rad/sec 

y L3 measured x distance from the GPIP normalized by u , sec 

y L4 MTjS measured z distance above the touchdown point normalized by U 0 , sec 

y L5 barometric altitude measured z distance above the touchdown point 

normalized by u , sec 

y L6 processed vertical velocity measurement normalized by U Q , non-dimensional 

y L7 processed acceleration measurement along the z-stability axis normalized 

by U Q , sec -1 

y L8 processed calibrated airspeed measurement normalized by u , non-dimensional 

y L9 acceleration measurement along the x-stability axis normalized by U Q , sec -1 

2, Lateral 

b 31 constant gain aileron servo position for lateral acceleration 

inovation, 1 /sec 

roll innovation, rad 

v 2 yaw innovation, rad 

v 5 y position innovation, sec 

lateral acceleration innovation, sec 

y 1 measured roll attitude, rad 

y 2 measured yaw from runway heading, rad 

y 5 measured y distance from runway centerline, sec 

y 6 processed lateral acceleration measurement, sec~^ 

X 1 estimated roll attitude, rad 

x 2 estimated yaw attitude from runway heading, rad 

x 3 estimated inertial sideslip, non-dimensional 

x 4 estimated roll rate about the x-stability axis, rad/sec 

x 5 estimated roll rate about the z-stability axis, rad/sec 
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estimated y distance from the runway centerline, sec 


x. . 
i-D 


w. 






W c 


w. 


-measured perturbed rudder position, rad 

the i th row j th column filter gain for aircraft state estimation 
estimated normalized gust component about the z— body axis, non-dimensional 
gust state variable, sec 1 

estimated rotational gust component about the z-body axis, rad/sec 
estimated rotational gust component about the x-body axis, rad/sec 
normalized steady state wind along the y-body axis, non-dimensional 
normalized shear wind along the y-body axis, sec 
the i th row j th column filter gain for wind state estimation 


x 3 


3. Euler Transformations 


EB 


EB ij 


J ES 


ES i j 


estimated roll attitude, rad 
estimated pitch attitude, rad 

estimated yaw attitue from runway heading, rad 
3 x 3 Euler transformation matrix, body to Earth axis 
the i^ row j ^ column component of the L matrix 
3 x 3 Euler transformation matrix, stability to Earth axis 
the i*'* 1 row j ^ column component of the L matrix 


4. General 


gain matrix for Kalman Filter aircraft states 
gain matrix for Kalman Filter wind states 


third-order complementary filter estimate of the x-coordinate in the MLS 
frame, ft 

third-order complementary filter estimate of the y-coordinate in the MLS 
frame, ft 

third-order complementary filter estimate of the z-coordinate in the MLS 
frame, ft 


7 



V e c constants for alpha-beta filters used to smooth aircraft rotational body 
rates, nondimensional y 


Guidance and Control Variables and Parameters 


l * _ Long it u d i na 1 


counter variable for the flare airspeed command easy-on 

incremental change per iteration of the counter variable for flare airspeed 
command easy-on p 

flare airspeed command easy-on value, non-dimensional 
perturbed airspeed command for flare, non-dimensional 
flare gain increase easy-on value, non-dimensional 

constant fiare gain increase factor on inertial alpha or normalized velocity 
along the z-stability axis 

constant flare gain increase factor on vertical position error 
constant flare gain increase factor on vertical position error 
flare x distance argument, sec 

tangent of the instantaneous flare flight path angle 

desired vertical acceleration along flare trajectory, sec -1 
derivative of h 1 ', sec - ^ 

desired perturbed altitude from the glideslope, sec 
desired perturbed pitch attitude, rad 

desired perturbed inertial speed alopg the x-stability axis. non-dimensional 
desired perturbed inertial alpha, non-dimensional 
desired perturbed pitchrate, rad/sec 

desired perturbed altitude from the nominal glideslope path (positive down), 

desired perturbed inertial velocity along runway centerline (positive in 
direction of landing), non-dimensional 

command pitch rate due to cross-coupling and second order effects, rad/sec 
C °second^order Effects, ^ec" 1 *» *> cross-coupling and 
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commanded acceleration along the z-stability axis, due to cross-couplinq and 
second order effects, sec" 1 

predicted velocity along runway centerline, non-dimensional 

commanded velocity along z-stability axis du^ to cross-coupling and second 
order effects and predicted velocity, sec 

commanded pitch rate due to trajectory changes, rad/sec 

commanded acceleration along x-stability axis due to trajectory changes, 
sec" 1 

commanded vertical acceleration due to trajectory changes, sec" 1 

-commanded pitch acceleration due to trajectory changes, sec" 1 

commanded vertical velocity (positive down) due to trajectory changes, 
non-dimensional 

glides lope easy-on value for vertical velocity command 
glideslope easy-on value for gain increase 
glideslope gain increase factor 
error from desired pitch attitude, rad 

error from desired inertial speed along x-stability axis, non-dimensional 
error from desired inertial angle of attack, rad 
error from desired pitchrate, rad/sec 

vertical position error from desired trajectory (positive down), sec 

error from trim thrust, K lbs 

error from trim throttle position, deg 

error from trim stabilizer position, deg 

2 

integral of the vertical error, e^ sec 

Lx6 


feedback 

gain 

for 

the 

. th 

i 

control 

and 

the 

. th 

error signal 

feedback 

gain 

for 

the 

. th 
i 

control 

and 

the 

. th 
3 

wind estimate component 

feedback 

gain 

for 

the 

. th 

l 

control 

and 

the 

. th 
3 

desired path state 

feedback 

gain 

for 

the 

. th 

l 

control 

and 

the 

. th 
3 

commanded rate error s 

feedback 

gain 

for 

the 

. th 
i 

control 

on 

the 

vertical path error 


integrator, z 

* INT 


0 


zPi 


cTjl 


cL2 


cL3 


The 


h Lx3 (i) 


h Lx4 * 1 * 




feedback gain for the i th control on the vertical gain programmed as a 
function of distance to touchdown 

elevator position command, rad 

stabilizer rate command, rad/sec 

throttle rate command, deg/sec 

elevator servo position command, deg 

throttle servo rate command, deg/sec 

desired altitude to engage flare mode, ft 

contains initial values of feedback gains, h T . at start of flare 
for i = 1 , 2, and 3 LX1J 

contains initial values of feedback gains, h r , 4 , at start of flare 
for i = 1, 2, and 3 * 

contains initial values of feedback gains, h T at start at flare 

for i = 1 , 2, and 3 xlb 


Av FR rate of change of desired airspeed during flare, ft/sec 2 

AV FD change in desired airspeed per controls update cycle during flare, ft/sec 

AV FL maximum allowed change to airspeed command during flare, ft/sec 

h 1/2P 1//2 of the desired x-coordinate distance traveled during the vertical 

acceleration period during flare, sec 

h FLR normalized flare height constant (HFLR F /U Q ), sec 

HFLRp selectable flare height trajectory constant, ft 

XFLR F selectable x-coordinate distance from glidepath intercept point (GPIP) 
to desired touchdown point, ft 


FLR normalized touchdown point distance (XFLRp/U Q ), sec 

T Y 0 tangent of the desired glideslope angle, non-dimensional 

T Y to constant equal to the tangent of the desired glideslope angle for touchdown 
of the final segment of the flare trajectory, non-dimensional 

Ax f a constant flare trajectory distance normalized by u Q (h 1 ^ 2p - x FLR ), sec 

x f flare trajectory distance parameter 

CWB n ,CWB 14 ,CWB 15 ,CWB 16 constants used to compute desired inertial speed deviation 

from estimated winds 
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x5 


~x6 


3 x7 


r INT 


INT 


INT 


x. . 
3-D 


C. . 

ID 


INTi 


INTi 


INT. 


u 


cl 


c2 

6 

ac 

5 

RC 

6 

RC 

h DC 

DECRAB 

A ez 

A<j> 

A^ 

a xi j 


LGN 


error in desired yaw rate about the z-stability axis, rad/sec 

error in the desired position deviation perpendicular to the runway 
centerline, sec 

rudder position deviation from trim, rad 

integral of the error in^desired position deviation perpendicular to the 
runway centerline, sec 

gain factor on error in desired position deviation perpendicular to runway 
centerline 

integral of the perturbed roll attitude, sec 

integral of the perturbed yaw attitude, sec 


feedback 

gain 

for 

the 

. th 

i 

control 

and 

. th 
D 

error signal 


feedback 

gain 

for 

the 

. th 

l 

control 

and 

. th 
D 

estimated wind 

component 

feedback 

gain 

for 

the 

. th 
i 

control 

and 

. th 
D 

command rate error state 

feedback 

gain 

for 

the 

. th 
i 

control 

on 

lateral position error integrator 

feedback 

gain 

for 

the 

. th 
i 

control 

on 

the 

roll integrator 


feedback 

gain 

for 

the 

i th 

control 

on 

the 

yaw integrator 



aileron position command, rad 

rudder rate command, rad/sec 

aileron servo position command, deg 

rudder servo rate command, deg/sec 

rudder servo position command, deg 

desired altitude to engage decrab mode, ft 

logical to indicate if decrab mode is engaged 

constant gain for decrab easy-on factor EZ1 

constant gain factor on desired roll-rate command 

constant gain factor on desired yaw-rate command 

constants gains to compute desired rate of change of the i th 

desired trajectory parameter, from the trajectory parameters, z 

logical to indicate if localize track mode is engaged 


PRECEDING page blank not 


niU-,:.- 


n 


d 31' d 4l' d 51 constant gains used in computation of desired lateral path trajectory 
parameters due to wind disturbances 

^F6 upper limit of localizer gain factor, G 

y k 

JPHINT logical to indiate if roll integrator is engaged 
^Lim maximum commanded bank angle during decrab, rad 
a a constant filter gain for aileron command 

u imin mainimum value of i lateral command 
u imax maximum value of i lateral command 

A1,A2,A3 hysteresis constants used in computation of aileron command 


C 


■'w 


w 


Hr 


SWTC 


at 


Ce ne_ra 1 

output aircraft state prediction matrix 
output wind state prediction matrix 
discrete control transition matrix 
feedback gain matrix for aircraft states 
feedback gain matrix for wind states 
feedback gain matrix for trajectory states 
feedback gain matrix for prediction states 
a ^ r craft state error vector 

one-half second easy-on factor used to turn on elevator or aileron 
at glides lope capture or localizer capture 

constant equal to the increment of time between discrete equation 
updates, sec 


command 


Prediction Variables 


1 « Longitu dinal 


Li 

4 > . . 
Lij 


T Lij 


Lwi j 


predicted value of the i th state component 


constants used to compute the i th predicted longitudinal aircraft state 
from the estimated aircraft state, 

<* transition natri* for 


predicted aircraft state update 

lenient of the i th row, j ^ colum 
predicted aircraft state update 


•lement of the i th row, ) th column of the wi nd transition ma t r i * for the 
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x.1 j 


W T ' 

Ll 


Lw. . 
- 13 


LI 


L5 


L6 


L7 


L9 


Lb. . 
- 13 

^L6 


f L7 


r L8 


7 L9 


"Lxij 


"Lwi j 


constants used to compute the predicted longitudinal aircraft state 
form the j cross-coupling and second-order command terms, n — 

predicted value of the i*^ 1 wind state component 

element of the i ^ row, j ^ column of the wind state transition matrix 

estimate of bias in the barometric altitude measurement, sec 

estimate of bias in the barometric altitude measurement, sec 

estimate of the bias in the vertical velocity measurement, non-dimensional 

estimate of the bias in the normal accelerometer measurement, sec” 1 

estimate of the bias in the longitudinal accelerometer measurement, sec 1 

element of the i^ row, j ^ column of the bias estimation matrix 

predicted vertical velocity measurement, non-dimensional 

predicted normal accelerometer measurement, sec 1 

predicted calibrated airspeed measurement, non-dimensional 

predicted longitudinal accelerometer measurement, sec” 1 

constants used to compute the predicted i th longitudinal measurement 
from the j predicted longitudinal aircraft state 

Vi 

constants used to compute the predicted longitudinal measurement 

from the j th predicted longitudinal wind state 


2, Lateral 

A 

predicted roll attitude, rad 

A 

X 2 predicted yaw deviation from runway heading, rad 

A 

X 3 predicted inertial sideslip, non-dimensional 

predicted roll rate, rad/sec 

/N 

x 5 predicted yaw rate, rad/sec 

predicted lateral position deviation from runway centerline, sec 
predicted rudder position from trim 

$ij element of the i th row, j column of the aircraft state transition matrix 

for the predicted aircraft state update 

Y j* element of the i row, j column of the control transition matrix for 

predicted aircraft state update 


< 


w. 


jo 


*. . 

id 


w. 


w. . 
- 1 D 
b . 


g * 

g v 

g f 

- y 


L6 


H 

A 

e 


'xi^ 


W1 j 


element of the i row, j th column of the wind transition matrix for 
predicted aircraft state update 

element of the i row, j column of the cross-coupling and second-order 
effects transition matrix for aircraft predicted state update 

. , . th _ . 

the 1 predicted wind state 

element of the i th row, j th column of the wind state transition matrix 

estimated bias in the roll attitude measurement, rad 

estimated bias in the yaw attitude measurement, rad 

estimated bias in the lateral accelerometer measurement, sec ^ 

gain for the roll bias estimate 

gain for the yaw bias estimate 

gain for the lateral accelerometer bias estimate 

predicted lateral position deviation from runway centerline, sec 

predicted rate of change altitude due to nonlinear effects, nondimensional 

estimated altitude above the touchdown point, ft 

predicted sideslip parameter, non-dimensional 

predicted sideslip parameter, non-dimensional 

command constants used to compute the predicted i th lateral measurement 
from the j predicted aircraft lateral state 

constants used to compute the i th predicted lateral measurements from 
the j lateral predicted wind state 


3. General 


ESij the * ^ element computed from predicted attitudes (<j>,6,t|0 of a matrix to 
transfer vectors from the stability axes to the Earth reference axes 

<(> predicted roll attitude, rad 

9 predicted pitch attitude, rad 

Y predicted yaw deviation from runway centerline, rad 


Initialization Parameters 
U Q selected reference calibrated airspeed, ft/sec 
WE I TO aircraft empty fuel weight, lbs/10 
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o 

Th 


HFLR C 


XFLR, 


YTD 


GSC 

a GSCi 

Pw- 


^Loc 
TFQ 
a 1 9 a 2 ' 


“k 


k 


a 


k 


selected qlideslope angle, rad 
trim angle of attack, rad 
trim pitch attitude, rad 
trim engine thrust, K lbs 
trim throttle position, deg 
trim stabilizer position, rad 

height of flare path above selected touchdown point, ft 

distance of selected touchdown point from the glideslope intercept point, 
tangent of desired flight path angle at touchdown 
calculated desired flare altitude, sec 

desired altitude to engage glideslope capture mode, sec" 1 

constants used in the computation of h 1 

GSC 

selectable constant having a value between 0 and 1 used for initializing 

predicted wind states w , and w ^ 

L4 L5 

desired distance from runway centerline to engage localizer capture 
mode, sec" 1 

total fuel quantity, lbs 

a^/b^i ,b 2 /b^, computed constants ued to determine trim inertial angle of 
attack, cc o , and trim thrust, T Q 

perturbed inertial velcoity along x-stability axis computed during 
pre-capture, non-dimensional 

computed inertial sideslope used during pre-capture 

computed inertial angle of attack used during pre-capture 


ft 


GENERAL DESCRIPTION OF SYSTEM AND TEST VEHICLE 

The DIALS is a complete digital autoland system consisting of a longitudinal and 
a lateral control law. The form of each control law is essentially the same and thus 
the functional blocks shown in figure 1 represent both control laws. Each control 
law consists of a fixed gain fullstate Kalman filter which processess innovations — 
the differences between sensor measurements and predicted measurements. Each con- 
tains a guidance algorithm which calculates guidance errors from the Kalman filter 
estimates and desired flight path trajectory commands (e.g., desired glideslope 
angle, desired calibrated airspeed). Finally, both control laws have a control 
algorithm which computes the servo commands needed to drive the guidance errors to 



zero. The commands of the longitudinal control algorithm are elevator position, 
ec' stabilizer rate, 6 gc , and throttle rate, <5 The commands of the 

lateral control algorithm are aileron position, <$ ac , and rudder rate, 6 R( _,. 

Each control law was designed independent of the other; but for the implementa- 
tion some cross coupling of the variables was used for each control law in the 
processing of the innovations and the desired flight path trajectory commands to 
account for some second order effects (e.g., the effect of roll on lift). 

Each control law has three control modes. Those of the longitudinal law are 
glideslope capture, glideslope track, and flare. These modes were designed to 
operate with standard and steeper glides lopes (the pilot selects a glideslope between 
2.5 and 5.5 degrees prior to system engagement). The control modes for the lateral 
law were localizer capture, localizer track, and decrab. The localizer capture mode 
was designed to provide satisfactory capture for ground intercept track anqles 
between 0 and 60 degrees of runway centerline. 


The test vehicle, shown in figure 2, was a small twin-jet commercial- type air- 
craft (ref. 8). it was especially equipped for flight test research of advanced 
navigation, guidance, and control laws and advanced displays. The aircraft was 
equipped with an aft research cockpit slightly forward of midship, triplicate flight 
control computers, a navigation computer, a display computer to drive electronic 
attitude and map displays, triplex inertial navigation systems and central air data 
computers, triplex sensors for the primary control variables, and a data recording 
system which stored in excess of 200 variables on magnetic tape. Of these variables 
forty-eight were software selectable internal flight control computer variables. 

Also located within the vehicle were test stations to monitor the computers, dis- 
plays, and recording systems during the flight experiment. 


CONTROL SYSTEM IMPLEMENTATION 


Sensor Measurements 


A total of 
The longitudinal 
used nine. Some 
were as follows: 


twenty-four sensor measurements were input to the control system, 
control system used seventeen of these measurements and the lateral 
signals were used in both systems. The signals input to the system 


Pitch (0), roll ( 4> ) , and yaw (i|> T ) attitudes were obtained from the Inertial 
Navigation System (INS) gimballed platforms. The pitch and roll were taken directly 

from the synchro outputs of the INS and thus were actually, respectively sin 0 
and sin <j> , 


Roll rate 
gyro’s. 


(p), pitch rate (q) , and yaw rate (r) were used from standard rate 


f z ) were obtained from a set 
The specific 


Specific forces in the aircraft body axes (f t , L , were c 
of 3 accelerometers mounted near the aircraft's center-of-gravity . 
force along the x body axis, ff x , was from the longtudinal accelerometer, ^and the 
specific force along the y body axis, f was from the lateral accelerometer. The 
output of the normal accelerometer has its output biased so that it has a zero output 
for steady level flight. In other words, the one g of gravity acting on the 
accelerometer in level flight was biased out. This accelerometer also gives a 
positive output for acceleration upward. Thus in order to get a right-handed 
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coordinate frame of specific forces for rotation to other axes and to properly 

account for gravity in the rotations, one g was added to the normal accelerometer 

output and then the negative of that quantity was obtained to give specific force 

along the z body axis, f . 

z 

Engine pressure ratios from the right and left engines, ERP1 and EPR2, were 
obtained from synchro outputs of engine sensors. 

Barometric altitude, hg, altitude rate, hg, and calibrated airspeed, v , were 
obtained from the central air data computers (CADC) • 

Radio altitude, h R , was input from the radio altimeter. 

The elevator servo position, <5 eg , was taken from the LVDT output of the 
elevator servo and the stabilizer position, s, was determined from the output of the 
trim potentiometer used for an autotrim algorithm of the test aircraft. 

The left and right throttle positions (Th Rp , Th Lp ) were from outputs of 
potentiometers mounted on the left and right throttle lever mechanisms. 

The right aileron servo position, , was taken from the LVDT output of the 

"A" hydraulic system and the left aileron §ervo position, , was taken from the 

LVDT output of the "B" hydraulic system. B 

Three MLS signals — the elevation signal, e£, the azimuth signal, Az, and the 
range or DME signal, R — were used as a position measurement. These signals were 
not input directly to the DIALS software. Instead, they were transformed to obtain 
x, y, and z components of a right-handed coordinate frame. These components were 
then input to the DIALS processing. The details of these transformations are given 
in appendix B. 

For reasons which will be discussed later, filtered vertical velocity from the 
MLS processing software, ^MLS was mac ^ e available to the pjALS software in place of 
hg by a software switch. The MLS vertical velocity is the output of a third-order 
complementary filter (see appendix B) which is driven by the z component of the run- 
way coordinate frame, mentioned above, and the vertical acceleration derived from the 
transformed body-mounted accelerometers. 


Discussion of Implemented Equations 

The DIALS system of equations (a modern multi-input, multi-output) control 
theory design) is a set of matrix equations plus the logical expressions required to 
engage and transition between the various control modes. For the computer simula- 
tions used during the DIALS design and development phase, these equations were imple- 
mented largely with matrix operation subroutines. However, for the flight test 
implementation, such an implementation was not possible due to flight computer limi- 
tations and real-time constraints. The flight computers were fixed-point 16-bit 
machines and were not suited to handle matrix operations efficiently. Thus, the 
matrix equations were broken down into individual element-by-element multiply and add 
operations. In addition, due to the fixed-point machine, variables were appropri- 
ately scaled to meet range and round-off error requirements — the method for determin- 
ing round-off error requirements is discussed in reference 7. 
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The specific equations implemented are given in appendix A. The appendix also 
includes the scaling of the variables, a flow chart of the equations and logic, and 
tables specifying values of the constants that were used in the test flights. The 
general equations which summarize the calculations done within the flight control 
computers are given below and are written in vector and matrix form for conciseness 
The calculations were executed once for the lateral control law and once for the 
longitudinal control law for each DIALS update cycle which was 0.1 seconds. The 
general equations are: 


Measurement Innovation Computations 


v - 
k 


- 


where y k is a processed measurement vector, y ts a predicted measurement vector, 
and v, is an innovation vector. * 


Steady-State Kalman Filter Computations 


\ * p v 


x k 


w. = w, + F v 
k k w k 


L ES ” l eb l bs 


where x 
vectors , 
vectors , 


and 


w. 


and 

and 


\ * are respectively the estimated and predicted aircraft state 
* are respectively the estimated and predicted wind state 

1 av*Ck rrnnn -ivt 4- * • . * 


vectors, F and F are constant matrices, L is the aircraft stability axis to 
earth axis Trunway frame) transformation matrix, E \ is the aircraft body axis to 
earth axis transformation matrix, and L the air£?aft stability axis to aircraft 
body axis transformation matrix. 


Guidance Variables Computations 


z k ' *z Vi + Vi 


C. = f(x,z) 


\ * f(x > 


e k x k “ z k 
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where z^ and C, are desired path command vectors , <f> is a matrix with time 

varying terms, is a vector containing non-linear second-order cross-coupling 

effects, and e^ is the trajectory error vector. 


Control Law Command Computations 


u. 


= -H e 
x k 


-H w 
w k 


H z v 

z k 


- H. 


( S - 


V 


+ added terms 


where is the control command vector; the added terms are integrator feedback and 

flare trajectory terms, and H , H , H and H ^ are matrices. 

x w z £ 

Prediction Computations 


Vi ■ % + r \ * r w ”k * *"’> 


w, „ = ♦ w, 
k+1 w k 


b = b + F v 
k+1 k b k 


Vi ■ C x Vi + C » Vi Vi 


where x^ and w^ are respectively the predicted aircraft state and wind vectors, 

b. is a vector of measurement bias estimates, and <b,r,r.¥,d),F,C, and C 
k w w b x w 

are matrices. 


All the DIALS equations were executed during two major frames within the flight 
computers — filtering and control commands in one frame and aircraft state and mea- 
surement predictions in the other. The DIALS was designed to be updated every 
100 msec whereas the time for one major flight computer frame was 49.152 msec or 
98.304 msec for two major frames. The effects of this time difference from the 
design time is small and is discussed in reference 7. 


Control Law Implementation Modifications 

Most of the modifications of the control law required for flight testing and 
those made to improve performance were accomplished in nonlinear batch simulations 
which were good representations of the test aircraft, engine dynamics, and control 
servos. However, one significant modification was made during ground testing of the 
flight software and another of less significance was made during flight testing. 
During real-time simulation ground testing in the Experimental Avionics System 
Integration Laboratory (EASILY), the elevator and servo commands were modified to 
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make the system insensitive to transport delays encountered in the EASILY Tests 
showed that some signals had an average delay of 100 msec with peak delays reaching 
200 msec. The delays arose because the digital flight control computers and digital 
aircraft simulation computer operated asynchronously and could not be synchronized. 
Although the maximum transport delays on the aircraft were only expected to be 
30 msec, the control laws were modified to make them insensitive to 100 msec 
delays. The insensitivity was achieved by a simple modification of the original 
elevator and aileron control commands while maintaining system stability (see ref. 7 
for analysis) as follows: 7 


<5 


c 


k 


a 



f 


u 

c 


where u c is the original unmodified command, is the previous modified 

command, 6^ the current modified command and a k_ i s a constant (a = 0.24758 for 

the elevator command and a = 0.170355 for the aileron command used for all test 
flights) . 

The modification made during flight testing was to replace the barometric sink 
rate signal (vertical velocity) from the central air data computer with an estimate 
of vertical velocity from a third-order complementary filter (see appendix B) . The 
modification was necessary because of the poor signal quality of the measured baro- 
metric sink rate for the test flights. The peak errors reached 30-40 ft/sec which 
was at least an order of magnitude higher than that used for design of the DIALS 
Kalman filter. The noise design values were more appropriate for cruise configura- 
lon of the test aircraft. The noise on this signal increases dramatically when the 
aircraft is m a landing configuration (flaps and gear down) and particularly in 
turbulent wind conditions. Thus to avoid redesign of the Kalman filter, the baromet- 
ric sink rate signal was simply replaced with the filtered vertical velocity signal 

after simulation tests showed satisfactory operation of DIALS with the filtered 
signal. 


During the EASILY tests, second-order digital prefilters described in appendix C 
were implemented m the flight computers to filter the body-mounted accelerometer 
measurements. These were implemented to prevent aliasing. In the early EASILY 
tests, aliasing was considered as a possible contributor to the instabilities which 

Jeft r in e th/?r d J° ^ T S6d the transport dela y s * However, the prefilters were 
left m the flight computers and were available upon activation of a software switch 

They were activated on one test run which will be discussed later. ' ' 

First-order low-pass software-selectable digital filters were also coded to 
remove noise from the body rate measurements during the EASILY tests. However the 
filters were not used in the flight tests since the EASILY tests showed that the 
filters caused a performance degradation after the modification for the transport 
delay was implemented. p 

^ rU ^° f ! fligHt tSSt ' several ^ ains were set ^ zero in the longitudinal 
control law. This test was made to verify simulation tests showing that zeroing the 

gains would not affect performance. The zeroing of the gains would result in elimi- 
nation of several multiply and add operations reducing code and real-time require- 
ments. These tests supported efforts to develop a "simplified" version of the DIALS. 
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FLIGHT TEST DESCRIPTION 


The flight tests were conducted at Wallops Flight Center (WFC) on Runway 22 (the 
runway which has an MLS installation)* Prior to a DIALS test run, the aircraft first 
proceeded to a calibration test point on Runway 22 (test point No. 12, 1000 ft from 
the stop end of Runway 22). The aircraft's vertical stabilizer, which had a laser 
tracking cube mounted on it, was positioned over this point to calibrate the WFC 
FPS-16 radar/laser tracking system. After calibration of the tracking system, the 
aircraft taxied to test point No. 74 for calibration and functional checks of the 
MLS. The aircraft's cockpit, which had the MLS C-band receiving antenna mounted 
above it, was positioned over this test point and the MLS signals were observed and 
recorded. The bias error on the MLS range or DME signal was computed and entered 
into the flight computer? then, on subsequent flight test runs, the bias was 
subtracted from the DME measurements. After the calibrations, the aircraft flew to 
the area to begin a DIALS test run. Upon arriving in the test area, the desired 
glideslope angle was entered into the flight computer (in a fully operational system 
the pilot would enter or select the desired glideslope). The pilot then positioned 
the aircraft at the desired flight level and ground track intercept angle specified 
for the test run using a velocity control-wheel-steering (VCWS) mode of the 
autopilot. (The general vertical and lateral path profiles for the flight tests are 
shown in fig. 3.) When the pilot achieved the desired flight level and ground track 
intercept angle, the pilot selected the desired reference calibrated airspeed for the 
landing on the mode control panel and, if not already engaged, engaged the calibrated 
airspeed hold mode. The desired ground track angle and flight level were then auto- 
matically maintained by the autopilot engaged in the VCWS made until engagement of 
the DIALS localizer and glideslope capture modes. After localizer and glideslope 
engagements, the aircraft flew to touchdown controlled automatically by the DIALS 
control laws. At touchdown, the system was automatically disengaged and the safety 
pilots then took control of the flight. The safety pilots then performed a nosewheel 
letdown and executed a touch-and-go in preparation for a subsequent test run. 


FLIGHT TEST RESULTS 

Selected time history plots of data from the test runs of the last five DIALS 
flights (R353, R355, R357, R361 , and R364) will be presented in this section to show 
the performance of each of the control modes of the DIALS--localizer capture, 
localizer track, decrab, glideslope capture, glideslope track, and flare. 


Localizer Capture 

The DIALS begins the capture of the localizer (runway centerline) according to 
its capture criteria logic. The criteria is a function of the distance of the air- 
craft from the centerline (Ay), aircraft true heading (Ai^), track angle (A^ TK )/ and 
estimated winds (w^). Localizer captures were performed in various wind conditions 
and track angles relative to the runway centerline. Figure 4 shows the capture 
performance for intercept angles from 20 degrees to 50 degrees which were performed 
on successive runs of test Flight R355. The data for deviation from runway center- 
line (Ay) shows that the capture overshoot increases with increasing intercept angles 
from no overshoot for a 20 degree intercept to about 50 feet for a 50 degree inter- 
cept. Figure 5 shows a 30 degree intercept angle capture approaching from the left 
side (facing runway) and figure 6 shows 30 degree intercept angle capture approaching 
from the right side, both in the presence of a 1 2 knot crosswind. The overshoot of 
the runway centerline is small in both cases and the yaw attitude (^) achieves the 
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proper crab angle in a smooth manner (without oscillation) which is attributed to the 
good estimation of the runway crosswinds . 


Localizer Tracking 

The localizer tracking mode maintains the aircraft on the runway centerline with 
the proper crab angle in crosswinds while holding the roll attitude at zero or wings 
level attitude. The plots of figure 5 illustrate the system's ability to maintain 
the above conditions and, in addition, show the response to a crosswind shear. The 
shear begins at 100 seconds into the run, and is characterized by the crosswind 
decreasing to zero over a small altitude distance (a shear wind of 8 knots/100 ft). 

A drift off runway centerline was induced by this shear of approximately 30 ft, but 
the aircraft was corrected back to centerline within 10 seconds after the shear- 
stopped. The good tracking performance of the system in the presence of gust winds 
is illustrated in the plots of figure 7. In the period of 60 to 95 seconds, the 
crosswinds magnitude changed four times — each change of approximately 10 knots. The 
maximum excursion from centerline during these gust encounters was 30 ft. 


Decrab 

The function of the decrab mode is to align the aircraft's true heading to that 
of the runway in the presence of crosswinds by performing a sideslip maneuver. 
Examination of the roll angle, <j>, and the heading or yaw angle relative to the 
runway, Aiji, (see figs. 4(b), 4(c), and 4(d)) near the end of the runs illustrates the 
maneuver and the variable, Ay, shows the good tracking performance achieved during 
the maneuvers for crosswind magnitudes of 6 knots. For high crosswinds of 18 to 
20 knots, the DIALS decrabed the aircraft, but limiting occurred on both the aircraft 
bank angle (limited to 5 degrees for safe clearance of engine pod and wing tip at 
touchdown) and the rudder position (see figs. 7(a) and 7(b)). The rudder limiting 
was set on the test vehicle as a function of dynamic pressure sensed in the elevator 
servo system to a value low enough to permit the safety pilots some remaining 
authority to override the automatic system commands. 


Glides lope Capture 

The objective of the glideslope capture mode was to capture the glideslope with 
low overshoot and low settling time for pilot-selected glideslopes between 2.5 to 
5.5 degrees. The last test flight (R364) contained a series of four consecutive runs 
in which the glideslope was increased from 3 to 5 degrees. The flight test data for 
3, 4, 4.5, and 5 degree glideslopes is shown, respectively, in figures 8, 9, 10, and 
11. The winds for these runs were high. The steady headwind component was' 10 to 
12 knots and the steady crosswind component was 18 to 20 knots. In addition to the 
steady wind, rapidly changing gusts of 8 to 1 0 knots were encountered (e.g., see the 
airspeed, V , of fig. 8(a) at 110 to 120 seconds). 

a 

The start of the capture is shown in the figures when the DIALS elevator 
command, u c j i -| ' anc ^ throttle rate command, u , become active. The capture 
completion occurs essentially when the glidesfope error. Ah, reaches zero for the 
first time. The capture time for all cases was less than 10 seconds and the maximum 
value of the first crossover of the glideslope error was 10 feet (see glideslope 
error, Ah, of figs. 8(b), 9(b), 10(b), and 11(b)). ‘ 
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In a non-gusty wind condition, the glideslope capture is normally initiated by a 

commanded reduction in throttle and a coordinated small positive commanded elevator 

deflection for pitchdown; however, in gusty conditions these control commands may 

change considerably due to the feedback of estimated winds and aircraft states to the 

control law* For example, note in figure 9, that there was no initial positive 

elevator deflection and that the throttle rate command became positive after a small 

initial negative value. In this case, the control law was adjusting for a decreasing 

headwind condition occurring during the capture (see V of fig. 9(a)). 

a 

The glideslope capture criteria was designed to adjust for different initial 
aircraft conditions. Figures 12 and 13 illustrate this adjustment capability. 

Figure 12 shows the capture of a 3 degree glideslope while the aircraft was descend- 
ing (approximately 4 ft/sec). Since the throttle was at the position required for 
the 3 degree glideslope tracking, the control did not command a throttle decrease. 
Only a small initial positive elevator command was required to initiate the capture 
which was completed smoothly without overshoot. Figure 13 shows a 4.5 degree glide- 
slope capture for the aircraft climbing (approximately 12 ft/sec) at capture initi- 
ation. For this case, a large throttle decrease was commanded along with a positive 
elevator command and, as before, the capture was completed smoothly with low 
overshoot. 


Glideslope Tracking 

The glideslope tracking mode maintains the aircraft along the selected glide- 
slope while controlling the calibrated airspeed to the chosen reference speed. The 
DIALS achieves airspeed control indirectly by commanding a desired inertial speed 
which is computed by subtracting the estimated headwind from the selected reference 
airspeed. The control is illustrated in the plots of calibrated airspeed, V , and 
desired inertial speed, U + z , of figure 11. At 73 seconds, a headwind gust 
occurs increasing the airspeed while, at the same time, the desired inertial speed 
begins to decrease. At 78 seconds, the throttle rate command responds to the 
decreasing desired inertial speed. Examination of the airspeed data in figures 8 
through 11 shows that its nominal value was 4 to 5 knots above the reference speed 
(the middle value on the y-axis for airspeed, V , is the desired reference airspeed) 
with peak deviations from the nominal of +5 knots. 

Glideslope tracking performance in wind gusts of 8 to 10 knots is shown in the 
plot of glideslope error, Ah, of figures 8 through 13. These plots can be compared 
to the plot of Ah in figure 14 which is typical of the glideslope performance for 
mild wind conditions. For the gusty wind conditions, the peak glideslope error devi- 
ations generally do not exceed 10 feet? but, in the data of figure 8, a peak glide- 
slope error of 20 feet occurred around 90 seconds apparently due to strong vertical 
up and down wind gusts that began around 90 seconds. Note that at 80 seconds, the 
aircraft is rising above the glideslope while both controls are in a direction to 
reduce the glideslope error. The elevator is going positive to pitch the aircraft 
down while the throttle is reducing. The airspeed, V , does not show a signiicant 
increase at that point, but the normal acceleration signal, a^, does show an upward 
acceleration indicating the presence of a vertical gust. Some sensitivity of the 
elevator command to certain gusts is shown in figure 8. At around 95 to 105 seconds, 
sharp headwind gusts were encountered (see V of fig. 8) which induced some 
oscillatory behaviour in the elevator command for 3 to 5 seconds. 

For mild wind conditions (relatively steady winds and low frequency gusts), 
represented by the plots of figure 14, peak glideslope errors do not exceed 5 feet 
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(see the plot Ah). The wind conditions for the data of figure 14 consisted of a 
6 knot headwind component and an initial 12-13 knot crosswind component which 
gradually reduced to zero at low altitude (350 to 400 feet). Note that, similar to 
the gusty wind conditions, the airspeed control for mild wind conditions also has 
about a 5 knot positive bias. 

The glides lope performance may have been affected by a lag in the throttle 
response to the throttle command. The lag which was identified in the flight data 
but not included in the design or simulation development occurred between the forward 
throttles, which directly drive the engine fuel flow valves, and the aft throttles 
which drive the forward throttles. The aft throttles are driven by the integrated 
throttle rate command, u . An examination of 6 of figure 14(c), which is 

EST 

the average position of the forward throttles, and the aft throttle position, 

°Th ' shoWS lags between the respective peaks and valleys of this data. For some 

peafcs and valleys, the lag can be characterized as a pure lag and for others it 
appears more like hysteresis. This lag can also be identified in figure 14(b) by 
examination of the throttle rate command, u cl3 , and the throttle position, 6-. . At 
120 seconds u CL3 is negative, but 6 Th does not move in the negative direction 
until at least a 1 -second delay. The lag, in general, could be characterized as a 
one second lag. The slow oscilatory motion of the throttle in figure 14(b) 

(“20 second time period) may be due to this unmodeled lag. Simulation studies which 
could quantify the effects of the lag have not been performed. 

The elevator activity before the DIALS glideslope capture (control by non-DIALS 
commands) and its activity immediately after capture appears to be about the same for 
the data of figures 8(b), 9(b), 10(b), 11(b), and 14(b) (see plot); but in the data 
of figures 12 and 13, more elevator activity occurs immediately after glideslope 
capture. Figures 12 and 13 represent data from an earlier flight (Flight R357, runs 
7R and 8R) than the data of the figures discussed above. On the earlier flight the 
gain on pitch rate feedback ( h Lx14 > was -2.0 compared to -2.2 for the later flights. 
The lower gain may have caused increased elevator activity; but the possibility of an 
unintentional increase in the glideslope tracking gains by an erroneous easy-on 
function seems more likely to have caused the increased activity. 


Flare 

The general objective of the flare mode is to transition the aircraft from the 
selected glideslope condition, whether a standard or steep one, to an acceptable 
touchdown condition. The DIALS flare mode achieves this condition by following a 
flare path from the glideslope to a specified touchdown point (a selectable param- 
eter) on the runway. Upon selection of the desired glideslope, flare path parameters 
are set to generate a desired vertical path profile (altitude versus distance) during 
the flare maneuver. During the maneuver, the flare control law commands the aircraft 
to pitch up to keep the nosewheel clear of the runway upon main gear touchdown; it 
also slowly commands a reduction in the airspeed rather than a reduction in the 
throttle as is done for existing flare laws. 

The performance of the flare law for mild wind conditions is illustrated in the 
data plots of figure 14 which shows a flare from a 4.5 degree glideslope. The plots 
of pitch, 0, and vertical velocity, h show a smooth flare maneuver which was 

engaged at 180 seconds. At touchdown, Sie pitch attitude is a positive 2 degrees and 
the vertical velocity is -2 ft/sec while the deviation from the flare path. Ah is 
nearly zero. The commanded reduction of airspeed during flare is achieved by 
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reducing the desired inertial speed, U q + z , and the resulting airspeed decrease 

can be seen in the plot of V . ° h 

a 

The performance of the flare law for gusty wind conditions is partially 
illustrated by the plots of figures 8, 9, 10, and 11; these are partial results since 
no landings were made for the runs of these figures because of lateral drifting due 
to high crosswinds and rudder limiting. However, prior to the go-around maneuver, 
the flare control, executing the flare from different glideslope angles, had 
established a positive pitch attitude and reducing vertical velocity for each glide- 
slope angle. A significant vertical deviation from the flare path is shown by the 
Ah plot of figure 8 which was induced in part by gusty conditions and in part by 
abnormal lateral control maneuvers due to the high cross winds and rudder limiting. 

Ten completely "hands-off M automatic landings were performed during the flight 
tests. For these landings, the mean touchdown point was 1164 feet from the glide- 
slope intercept point compared to a commanded value of 1300 feet (see ref. 7). The 
standard deviation from the mean touchdown point was 244 feet. The mean touchdown 
vertical velocity was 2.4 ft/sec with a standard deviation of 0.74 ft/sec. The 
commanded touchdown vertical velocity is a function of the ground speed which for a 
125 knot ground speed is -2.2 ft/sec. 


CONCLUDING REMARKS 

The implementation of a modern control theory, direct-digital designed, full- 
state feedback autoland system on a 16-bit fixed-point flight computer was achieved 
and was successfully flight tested. The DIALS thus became the first modern control 
theory autoland design to be successfully flight tested on a small commercial type 
jet aircraft. The tests demonstrated that the methods used are a viable approach to 
the design of multiple-input multiple-output systems which can provide increased 
performance through the coordination of the multiple control commands. The system 
also demonstrated expanded capabilities over current systems. The expanded 
capabilities include the capture of glides lopes and localizers with lower overshoot 
and settling time, the capture of steeper glideslopes, and the ability to perform the 
flare maneuver from various glideslopes. 

The DIALS achieved low overshoot captures (average value of 24 ft with standard 
deviation of 26 feet) of the runway centerline for capture angles of 20 degrees 
through 50 degrees. The aircraft was generally settled and tracking the runway 
centerline less than 40 seconds after capture initiation which is a factor of 3 less 
than current autoland systems. 

The DIALS maintained small glideslope and localizer tracking errors for mild 
wind conditions — steady winds of 10 knots and gusts of 3 to 5 knots. For these 
winds, the glideslope and localizer errors did not exceed, respectively, ±5 feet and 
±20 feet. 

For gusty wind conditions, steady winds 20 to 25 knots and gusts 8 to 10 knots, 
glideslope and localizer tracking errors generally did not exceed, respectively, 

±10 feet and ±30 feet. However, on some occasions, the longitudinal control law 
exhibited an oversensitivity to gust wind conditions resulting in oscillatory 
behavior. Pitch changes of 8 degrees within 5 seconds and glideslope errors of 
20 feet occurred. Of course, the aircraft must make pitch changes, in the absence of 
direct-lift control, to tightly track the glideslope in gusty wind conditions. The 
pitch changes that did occur were in the right direction to correct the glideslope 
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errors, but appeared to be too large or, in other words, an over correction. Another 
factor that contributed to the pitch behavior in these gust conditions was a signifi- 
cant lag of approximately one second in the response of the throttle to a command 
which was not accounted for in the system design nor in non-linear simulation 
development. 

The DIALS performed successful decrab maneuvers in crosswinds up to 6 knots 
Decrab maneuvers in crosswinds of 18 to 20 knots were performed but lateral drifting 
occurred due to rudder and bank angle limiting. (Rudder limiting occured at approxi- 
mately 10 degrees for approach condition on the research system to permit safety 

pilot override.) Crosswinds between 6 and 18 knots were not encountered in the 
tests. 


The longitudinal capture criteria and capture control mode demonstrated an 
insensitivity to initial conditions and the selected glideslope angle. Various 
glides lopes-- 3 degrees to 5 degrees— were successfully captured with low overshoot 
(typicaily less than 5 feet) without system gain changes from descending, level, and 
climbing initial flight conditions. 


The flare mode successfully flared the aircraft from 3, 4, 4.5, and 5 degree 
glides lopes with automatic landings resulting for 3 and 4.5 degree glideslopes 
Successful landings for the other glideslope angles would have been achieved (pitch 
attitude and vertical touchdown velocity were acceptable), but due to lateral drift 
in high crosswinds, go-around maneuvers were executed. 

In general, the flight tests demonstrated that each of the six control modes 
successfully performed its intended function. Improvements to the performance of 
each could most likely have been achieved given additional flight tests and data 
analysis time. For instance, no tests were performed for different flare trajec- 
tories; that is, different selectable trajectory parameters in the flare control law 
such as the desired touchdown point were not tested. No time was spent to fine tune 
these parameters to the aircraft from values set during the non-linear simulation 
development. 


The implementation and flight tests of the DIALS illustrated the importance of 
accounting for system transport delays and lags. The transport delays encountered 
during the EASILY ground tests reguired a modification to the DIALS before the ground 
tests could be successfully completed. The modification may not have been required 
or the flight tests since the average transport delay on the aircraft was estimated 
to be a factor of 3 less than that measured in the EASILY tests. Thus, the message 
is that the ground test facility should have transport delays which are on the order 
of that of the flight system. The lags encountered between the aft throttle position 
and the forward throttle position, which were on the order of 1 second, may have 
contributed to some throttle instability, especially in gusty wind conditions. The 
message m this case which is not new, but which was illustrated in the data, is to 
make sure that all significant lags are accounted for in the design and particularly 
m simulation testing and development. 9 
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APPENDIX A 


DIALS SOFTWARE EQUATIONS 7 SYSTEM FLOWCHART 
1 • Software Equations 

The following equations were implemented in the flight control computers of the 
test aircraft. For each equation and above each variable in the equation, there is 
the letter "B" followed by an integer number which specifies the scaling of the 
variable in the computer. For example, B-2 means that the maximum value of the 
variable is 2 or 0.25 or B6 means that the maximum value of the variable is 2 6 or 
64. From this scaling specification, the resolution of the variable can be found 
given that the computer was a 16-bit machine (one bit for sign). For the above 
examples the resolutions are, respectively, 0.25/32768 or 7.6294 * 1 0' 6 and 64/32768 

°*\ 1,9531 X . 10 (2 = 32768) ’ for some variables double precision was used and 

this is indicated by DP enclosed in parentheses following the scaling specification 

(e.g., B1 8 (DP). For this example, the maximum variable value is 2 18 or 262144 and 
the resolution is 2 18 /2 31 or 1.2207 * 10" 4 . 

pi- J hS values of ^ gains and constants used for runs 4, 6, and 7 of the last test 
flight are given in the last section of this appendix. 
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B . Filtered Longitudinal State Vector 
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C. Filtered Longitudinal Wind Vector 
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D. Lateral Innovations 
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B- 3 BO B3 B-8 


B-2 (DP) 
EZ5 


h < . . 
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Notes : 


B1 


B-2 


B1 B-2 
'o Z L 3. 


BO 


B1 


^ T( - = Z T1 ~ SCY + (T - T ) X 

Lb, LI, o L3. y a yd 


1 • G =0 for all flights 
4 

2. CWB15 should have been CWB14; the values for CWB14 & CWB1 5 are 
respectively, 0.99758 and 0.99452 (software specification error) 


B-2 

B-1 

BO 


x 

(cos ( 

B-2 

B- 

2 

?L1 k = 

" H 

B-2 

B-2 

B-1 

° L 2 k = ‘ 

x 

x , + 
L4 k 

B-2 

B-2 

B-2 

H = 


•x 

B-2 

B-1 

B-1 

X* 

x 

X' 

B-2 

B-2 


" L3 k = 

" nL3 k 



BO 


BO B-1 
: x_ 


BO B-1 


X 3 \ 

k k 


B-1 BO 
sin $ 
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B1 


BO 

BO 

BO 

B1 

X 

= (1 

+ T 1 

ES3 1 Y 

o 

ES1 1 5 

(1 + 



BO 

BO 

BO 

B-2 


+ 

(1 ES33 + T Y 

o 

1 ESI 3* 

X L3 



BO 

BO 

BO 

BO 


+ 

^ES32 + T Y 

o 

1 ES12 ) 

\ 


B1 B1 B1 

C L6, = C L6, ~ °L6, 
k k k 


B-1 



k 


B-2 
+ X L1 


k 


B-2 



B1 

SCY 

o 


If EZ7 < 1 


B1 ( DP ) B1 ( DP) B1 ( DP ) 
EZ7 = EZ7 4 * AEZ7 


and 


B1 (DP) 



B1 (DP) 



B1 (DP) 
* EZ7 


8. Glideslope Track Engage Logic 


If j gsC = 1 ' then Go to 9, 


If 


B8 B2 B1 0 


U 

o 


'Lx6, 


> 15 , 


then Go to 10* 


J GSC = 1 which means glideslope track engage 
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9. Easy-On Glides lope Gains 


B1 


If: j c _ 1 * 1, then Go to 10. 


B1 B1 (DP) B1 (DP) 

EZ4 = EZ4, + AEZ4 
k k-1 


B2 B2 B1 


B1 


DFGT = [l + (GF3 - 1) EZ4 ] 

k 


h . = h 

Lxi3 Lxoi3 


h r = h 

Lxi6 Lxoi6 


B2 

DFGT 

B2 

DFGT 


A 


> for i = 1,2,3 


10. Trajectory Errors 


B-2 


B-2 B-2 


e Lx1 X L1 Z L1 

k k 


B-1 
e 


B-1 B-1 


Lx2 k = X L2 k - Z L2 k 


B-2 
e 


B-2 B-2 


L * 3 k = Xl \ ' h 


B-1 

e 


B-1 B-1 


Lx \ = H ' H 


B2 

e 


B6 (DP) BO B10(DP) B6(DP) 

z 


Lx6 X L6 + T Y X L5. 
k k o k 


L6, 
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I 



B5 B5 

e = x „ 

Lx7 'Ll 

k k 

B6 B6 

6 Lx 8 = X L8, 
k k 


B-2 B-2 

e Lx9 = X L9 
k k 


1 1 . Altitude Error Integrator Logic 
If K FL = 1 skip to 12. 


B2 


If 



B8 

U 

o 


B1 0 

c z T . 

Lim 


then K__ 
FL 


1 


Else 


KFLC = KFLC - 1 and if KFLC < 0 then K pL = 1 
Else Go to 13. 


12. Altitude Error Integrator 


B4 ( DP) B4 ( DP) B2 ( DP) B-3 


z = z + e At 
INT INT Lx6 
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13. Flare Touchdown Parameters 

if J Flare = ^ en Go to 14. 

B5 BIO (DP) BIO (DP) 

X Lim X FLR ~ X L5 

k 


B5 


B5 


If X T - < °* 5 ' then x r • =0.5 

Lim Lim 


BO 

FL 


ELX 


B2 ( DP ) 

(2) e 

Lx^ 


B5 

x . 
Lim 


14. Longitudinal Commands 


B7 B9 

U CL1 = _h Lx1 ' e 
k 


B-2 

B8 

B-1 

B9 

B-2 



Lxl 

k 

-h 

Lxl 2 

e Lx2, 

k 

^Lxl 3 

e Lx3 

k 



B-1 

B5 

B2 

B2 

B5 

B1 

B6 

Lx4, 

k 

-h 

Lxl 6 

e Lx6, 

k 

” h Lx1 7 

®Lx7 ' 
k 

-h 

Lxl 8 

e 

Lx8, 


B9 

B-2 

B9 

B-2 

B9 

B-2 

B9 

B-2 

-h 

Lxl 9 

e 

Lx9, 

k 

-h 

Lwl 1 

w 

L1 k 

-h w 

Lwl 2 L2. 

k 

-h w 

Lwl 3 L3_ 

k 

B8 

B-1 

B8 

B-1 

B9 

B-2 

B9 

B-2 

~ h Lw1 4 

w „ 
L \ 

^Lwl 5 

w ^ 
L \ 

-h 

Lwl 6 

w _ 
L6 k 

-h 

Lwl 7 

w „ 
L \ 

B9 

B-2 

B8 

B-1 

B9 

B-2 

B8 

B-1 B5 

_h Lz1 1 


-h 

Lzl 2 

" L2 k ' 

’ h Lz1 3 

H ‘ 

^Lzl 4 

z — h 

L4 Lzl 6 

k 
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B9 

B-2 

B9 

B-2 

B9 

B-2 


“ h LC11 C L1 

k 

” h LC1 2 


” h LC1 3 



BIO 

B-3 

B6 

B1 

B3 

B4 

B7 BO 

“ h L^14 5 L4. 

k 

~ h LC16 


-h 

ZT1 

z 

INT 

-f h FL, 

c zPI 


B7 


u . . for ii . < u , „ 

Lminl CL1 Lmxnl 

k 


U CL1^= U CL1^ ^° r ^niinl* U CL1 ^ 'Ymaxl 


u T , for 11 . > u 

Lmaxl CL1 Lmaxl 

k 


B2 BO B2 


B2 


U = CX u + U 

CL1 e CL1 CL1. 

k k-1 k 


B2 


B4 B- 2 B3 B— 1 B4 B-2 B3 B-1 


U CL2 k h Lx21 S Lx1 k “ h Lx22 ®Lx2 k \x23 ®Lx3 k ~ h Lx24 ®Lx4 k 


BO B2 B-3 B5 B-4 B6 B4 B-2 

Lx26 S Lx6 ^Lx27 G Lx7 ^Lx28 S Lx8 ^Lx29 G Lx9, 

k k k k 


B4 B-2 B4 B-2 B4 B-2 B3 B-1 

•h W —Vi w — H Xfj —K t fj 

Lw21 Li. Lw22 L2, Lw23 L3, Lw24 L4, 
k k k k 


B3 B-1 B4 B-2 B4 B-2 B4 B-2 B3 B-1 

' h Lw25 W L5 " h Lw26 W L6. _h Lw27 W L7. " h Lz21 Z L1 _h Lz22 Z L2 

k k k k k 


B4 B-2 

B3 B-1 

BO B2 

B4 

B-2 

h _ z 

-h z 

-h z 

-h 

r 

Lz23 L3, 
k 

Lz24 L4_ 
k 

Lz26 L6, 
k 

LC21 

LI 
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B4 B-2 B4 B-2 B2 B-3 

~ h LC22 C L2 "V23 C L3, _h LC24 C L4, 
* k k 


B1 B1 B-2 B4 


B2 BO 


LC26 ^zT2 Z INT ^c h zP2 FL ELX 


r 


B-4 

u 

CL2, 


Lmin2 

for 

U CL2 < 
k 

U Lmin2 

l CL2. 

k 

for 

U Lmin2 

< u 

CL2. 

Lmax2 

for 

"C12 > 

k 

u 

Lmax2 


SWTC ( u cli + GSTAB Uq^) 57. 3/. 86) (deg) 
c 

where SWTC goes from 0 to 1 in one-half second (Easy-on) at glideslope capture. 


CL 3, 
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BIO 

B-2 

B9 

B-1 

BIO 

B-2 B9 B-1 

-h 

Lx31 

e 

Lxl , 
k 

-h e 

Lx32 Lx2 

k h Lx33 S Lx3 k h Lx34 S Lx4 

B6 

B2 

B3 

B5 

B2 

B6 BIO B-2 

~ h Lx36 

e 

Lx6, 

k 

h Lx37 S Lx7 

" h Lx38 e Lx8, \x39 ®Lx9 
K K 

BIO 

B-2 

BIO 

B-2 

B1 0 

B-2 B9 B-1 

-h 

Lw31 

X' 

’ h Lw32 

x 

^Lw33 

W L3. - h Lw34 W L4 
k k 

B9 

B-1 

BIO 

B-2 

BIO 

B-2 BIO B-2 

^Lw35 

w _ 
L5 k 

^Lw36 

w _ 
L6 k 

* h Lw37 

W tX ~ h Lz31 Z L1 
k k 

B9 

B-1 

BIO 

B-2 

B9 

B-1 B6 B2 

_h Lz32 

X" 

^Lz33 

X 

” h Lz34 

z — h z 

I>4 Lz36 L6 

K k 

B1 0 

B-2 

BIO 

B-2 

B1 0 

B-2 

~ h LS31 

x- 

’ h L<;32 

' L2 k ' 

" h LC33 



I 



B7 

U CL3, 


m i 

B-3 

B7 

B1 

B4 

B4 

~ h LC34 

' L4 k 

’ h LC36 

S - h 

zT3 

Z INT 

r 






U . _ 

for 

u < 

u 



Lim3 

CL 3 

Lim3 



/ U 

for 

u . < 

u 

< 

u 

1 CL3 k 


Lim3 

CL3 

k 


Lmax ■ 


B2 
f h 


BO 

FL 


u * for 11 > u 

Lmax3 CL3. Lmax3 

k 


B7 B7 

^Thc = U CL3 = t ^ irottle ra te command, deg/sec 
k k 


B. Lateral Control 


1 . Desired Path 


If DECRAB — 1, then Go to 2. Estimated Sideslip Parameters 


BO BO BO 


BO BO 


B1 B-1 B1 B-1 




B 

-1 

B1 

B1 

B-1 

B-2 

B2 

BO BO 

+ 

1 1 

z1 k-i 

+ ip 
r 1 3 

?z \-i 

+ *14 

Cz4 k-1 

4>. f - C c 
15 25 k-1 

B2 B2 

BO 


B2 

B2 

BO 

B3 

B-1 

Z 2 = *821 Z 1 
k 

+ 

k-1 

V, 

+ <J> 

z23 

Z3 k-i 

+ <f> 

z24 

V, 


B3 

B-1 

B1 

B1 

B3 

B-1 

B3 B-1 

+ 

d> 

z25 


+ *21 

Cz \-i 

+ ip 
Y 22 


+ U» r 

23 z3 
2J Z3 k _ 1 


BO 

B2 

B2 

BO 




+ 

ip 

24 

V 

k-1 

+ *25 

? z5 

k-1 
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B- 1 B- 1 BO B- 1 BO 


\ ' * 241 * * 2 « * 


4. r z4l ‘I. . ' r z43 ‘3. . ” *z44 z 4. . * *z45 z 5 


BO B-1 BO B-1 

k-1 'k-1 


B-2 B1 


BO B-1 B-3 B2 


B-1 BO 


^ ?z1 k -1 + * 43 ?23 k _-, + * 44 C *4 k _ 1 + *45 C z5 k _ l 


B-1 B-1 BO B-1 BO BO B-1 BO B-1 

z. 


•5 ’,51 V . * *z53 2 3 , * *z54 2 4 v , + *z55 z 5„ 

K ~ 1 k- 1 k-1 


k-1 


B-2 B1 
+ C 


BO B-1 B-3 B2 B-1 BO 


51 '•zl, . ^53 *"z3, . + ^54 ^z4, . + ^55 ^z5 


k-1 


k-1 


k-1 


k-1 


B6 B6 BO B4 B2 B6 BO 

2 6,_ ■ *461 2 1. . * *z62 2 2. * *■ 


B7 B-1 


'k — *k-1 262 "V, ' ’463 + *z64 Z 4 k _, 


B7 B-1 B6 B-1 B1 B1 B-1 B1 B-1 

>265 \-, + V/ * 62 H-, + *« S z3 


k-1 


B-2 B2 BO BO B-2 B2 

) r 
66 ^z6 


+ ,| ' 64 5z4 k _ 1 * *65 5 z5 k _ t 


* *66 5,6 )/ 2 6 


k-1 


2. Estimated Sideslip Parameters 


B-1 B-1 BO 


S k W 5 ^ EB22 
k 


B-1 B-1 B-1 B-1 

\ = \ ~ V “ V 

k k 


BO B-1 

Z 3 = _ " 6 k 
k K 
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If DECRAB 


1 


Go to 4 


B-1 

BO 

B-2 

• 

%■ 

_1 EB22 

\ 


3. Decrab Engage Logic 


If h DC 001 + X L6 < °' then ski P to (Lateral Path Command) 

k 


Else 


DECRAB = 1 , 


and 


EZ1 = 0, 




4. Decrab Roll Limit 


BO BO 


If > 


Lim 


Go to 6, otherwise. 


5. Decrab Trajectory 


BO BO B-1 BO 

Z 3 = Z 3 U + K EZ1 
k k 


B1 


BO B1 


B1 


EZ1 , = A EZ1 , „ + AeZI where 

k ez k-1 


If EZ1, > 1, then EZ1 V = 1. 


B1 BO B2 B1 
*1k ' A <|> *2k EZ1 k 


B-1 B-2 B2 B1 

*2k = \ X 2 ™k 
k 


6. Lateral Path Command 


B1 

C 


B2 B-1 B2 B-1 B1 




B-1 BO B-1 B1 B-1 B-1 

£ o = _a z „ - a „ z„ t z 

2 \ 1124 \ 225 5 k \ 


B-1 


z3. 


B- 1 BO B-1 BO BO B- 1 BO 


= -a 


B-1 B-1 


x31 1 ' x33 Z 3 V " a x34 Z 4, “ a X 35 Z 5, + Z - 


BO BO BO B1 B-1 B1 B-1 

^ z4 k a x43 Z 3 k " a x44 Z 4 k “ a X 45 Z 5 k 


BO BO BO B1 B-1 B1 B-1 

C c = -a __ z„ - a 


z5 


x53 3 x54 4. x55 5, 

K k k 


B2 B2 BO BO B2 B2 BO 


z6. 


= -a„,. z. a 


x61 1 x62 2, x63 *3, 


B1 BO B-1 

C<j>sS Sg 

\ vir + ^/ r * 

O 


B1 B1 B1 


c 


ul 


k 



+ C , 


k 


S6 


i 



B-1 


BO . B-1 


c<|> 


V --T +sc e ] r k 

Co O 


B-1 B-1 B-1 


x ■ + s 


B-1 BO B-1 B-1 B-1 B-1 B-2 

\ ’ ^31 ? k * X 5„ X L2„ ' x 4„ X L3 
k k k k k 


B-1 B-1 B-1 

C , = C - + C. 

u3 z3 3. 

k k k 


B2 B3 B-1 

5 4 = ~ d 41 ^k 

k 


B2 B2 


B2 


Q A = C „ + C. 

u4, z4, 4, 

k k k 


BO B1 B- 1 

\ ■ - d 5, 5 k 

k 


BO BO 


BO 


C _ - t, c + 

u5 k z5 k 5 k 


B2 


BO 


BO BO 


B- 1 B2 BO 


C = -(1 -1 ) x, - 1 (1 + X. „ ) + x^ CY 


'6, ’ ES22 3. ES21 

k k 


L2. 2, o 

k k 


BO B-2 BO BO 


-1 X X S 

ES23 L3 1, a 
k k o 
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B2 

B2 

B2 


^ U6 k = 

C z6 

k + \ 

7. 

Tra j ectory 

Errors 


BO 

BO 

BO 


e = 

x1 k 

\ 



B2 

B2 

B2 



\ 



BO 

BO 

BO 


e -v 

V 



B-1 

B-1 

B-1 


e x4 

X \ 

\ 



B-1 

B-1 

B-1 


^k = 

\ 



B6 ( DP) 

B6 ( DP ) B6 ( DP ) 


e = 
X k 

\ 



B-1 

B-1 




\ 


8. 

Localizer Track 

Engage Logic 

If 

J LGN =1 ' th6n 

skip to 9. 

If 

58 

( |« 6 I * u o < 

f G =1.0 

1^ J LGN 1 



Otherwise C LGN = C LGN - 1 ; If C LGJsJ > 0, Go to 12. 



9. Cross-Track Error Integrator 


B6 ( DP ) B6 ( DP ) B6 { DP) B-3 

y = y + e At 

J INT J INT „ x6, 

k k-1 k 


10. Easy-On Cross-Track Error Gain 


B3 B3 B3 B3 B3 

G = G + 0.05 if G < h 


y k y k-i 


Vi F6 


Else = h p6 


B6 ( DP ) B3 B3 


e = G e 
x6 y. x6 
k k k 


11. Roll and Yaw Integrators 


If (JPHINT = 0 and DECRAB = 0) then 


If | x 57.3 | < 1.0 deg then JPHINT = 1 


(JPHINT = 1 

and DECRAB = 

0) then 

B2 ( DP ) 

B2(DP) 

B0 

B-3 

<t> = 


+ 

At) 

INT 

k 

INT k-1 

Ik 
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If DECRAB =1 then 


B2(DP) B2(DP) B2 B-3 
^ INT k ~ ' , ' INT k _ 1 + X 2k AT 


12. Lateral Control Commands 


B7 ( DP ) B7 BO B5 B2 B7 BO B8 B-1 


B8 B-1 B1 B6 B8 B-1 B8 B-1 

■h e -h e -h e -h w 
X 15 X5 k X 1 6 X6 k X 17 X \ *1, \ 


B9 B-2 B9 B-2 B9 B-2 B8 B-1 

' h W W 2 _h w W 7 ~ h „ “ h w c 

12 \ 13 \ W 14 \ W 15 \ 


B9 B-2 B6 B1 B8 B-1 B8 B-1 
-h w -h (, -h C -h C 

”16 \ Si ',2 " 2 k S3 


B5 B2 B7 BO B5 B2 B1 


B6 B5 


C 14 U \ C 15 ?u5 k ^ U6 k hy iNT1 yiNT \ NT 1 ^ INT 


B5 B2 

-h 

INT 

INTI 


. for u < u 
Imin cl 1 

k min 

u for u < u < u 

c \ Imin cl u 1max 


u for u > u 

lmax cl, Imax 

k 


u = u + a u 

Cl k » <S-1 


TEMP = u „ 180+6 

cl. a 

k tt trim 



t> = SWTC 
ac 


TKMP_ 

~ B1 R2 

I TEMP I* A1 + A2 


If TEMP < A3 


If TEMP > A3 


where SWTC goes from 0 to 1 in one-half second (Easy— on) at localizer capture. 


B5 ( DP) 

B5 

BO 

B3 

B2 

B5 

BO 

B6 

B-1 


-h 

X 21 

6x1 k 

-h 

22 

6x2 k 

-h 

X 23 

“ x3 k 

-h 

X 24 

S x4 


B6 

B-1 

B-1 

B6 

B6 

B-1 

B6 

B-1 


-h 

X 25 

e *\ 

-h 

X 26 


-h 

X 27 

6x7 k 

-h 

21 

\ 


B7 B-2 B7 B-2 B7 B-2 B6 B-1 

-h w -h w -h w, -h w r 
w 2 w 3 w 4. w^_ 5, 

22 k 23 k 24 k 25 k 


B7 

B-2 

B4 

B1 B6 

B-1 

B6 

B-1 

-h 

W 26 

\ 

~ h c 

S 21 

C , -h 
U, k 4 22 


z3 


B3 

B2 

B5 

BO B3 

B2 

B-1 

B6 

24 

^ U4 k 

^25 

5 , -h 
U5 k 4 26 

^ U6 k 

-h y 

y INT 

y INT2 

B3 

B2 B3 

B2 





6 -h il> 

> INT & INT 

INT2 INT2 



for 

u < 

u 

2min 


c2 k 

2min 


for 

u . < 

U 1 

c2 k 


2mm 

c2 k 


for 

u > 

u 

2max 


c2 k 

2max 


B15 B2 

6 = u 

R °k ° 2 k 


B1 3 
180 
7T 


(deg/sec) 



(deg) 


6 = f 6 +6 

Rc J Rc, R 

k trim 


Note: The 6^ signal was integrated at the flight control computer 


major frame 


rate ("20 cps) 


TV. PREDICTIONS 


A. Longitudinal Predictions 


1 . State 


B1 BO 


B-1 B-1 BO B-2 BO 


\5 k 1 ESI 1 0 + X L2 r ) " X L2 k ° Y o + X L1 k SY o 


BO BO B-2 BO BO 

ESI 3 ' S V X L3* ^ ESI 2 X 3 

k k 


BO 

n 


BO 


B-1 


B-1 BO B-2 BO 


L6 k - ’eS3, <’ * \ 2|t » * * - L , k cr o 


B-2 

A 

*L1 


k+1 


BO 

BO 

B-2 

BO BO 


(1 ES33 

- cy ) 

o 

H 

+ 1 X 

ES32 3, 
k 


B-2 

B-1 B-1 


BO B-2 

B-1 

H + 

^L12 X L2 

k 

*L13 *1.3,. * 
k 

<P 

LI 4 

B-7 

B5 

B-8 

B6 BO 

B-2 

+ <f> 

LI 7 


* 

LI 8 

X + <±> 

L8 Li 9 

k 


B-4 

B2 

B2 

B-4 



B-1 


+ T u + y n 

L11 CLI^ Li 2 U CL2 k 
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> to 


BO 

B-2 

BO 

B-2 

B-1 

B-1 

S, 

“'•V 

+ Y r 

LW 13 

•«* 

+ Y r 

LW 14 

w 

M k 

B-1 

B-1 

BO 

B-2 

BO 

B-2 BO B-2 

YLW 15 


+ 

LW 16 

X 

+ i|> 
r L1 1 

n + Ui n 

LI, r L1 2 L2, 

k k 


BO 


B-2 


+ ^ n 

LI 3 L3, 


B-1 

A 

B1 

B-2 BO 

B-1 B1 

B-2 BO 

B-1 

X L2 

k+1 

= 4> 

L21 

*L1 + *1.22 

k 

X L2 + ^L23 
k 

X L3 + ^L24 
k 

*L4. 


+ 


B-6 B5 


L27 



B-7 

B8 

B1 

B-2 

+ <j> 
y L28 

x ^ + 

L8 k 

<b 

r L29 

X L9. 


+ 


B-5 


L21 


B2 B-1 0 B7 B1 B-2 

U CL1* Y L23 U CL3 + Y Lw W L1 

k k 21 k 


B1 
+ Y 


B-2 


BO 


w ^ + Y 

“2 3 L V “24 


B-1 

BO 

B-1 

w . 

+ Y 

w _ 

L4 k 

LW 25 

L5 


B1 
+ Y 


B-2 


B1 


B-2 


B1 


B-2 


L2 W L6 + *L21 \l + ^L22 

26 k k k 


B1 B-2 


+ 


L23 
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B-2 

* 

BO 

B-2 

B-1 

B-1 

BO 

B-2 

B-1 

B-1 

X L3 = 

k+1 

: <p 
y L31 

x + 

<p 

L32 

x t 2 + 
L2 k 

<p 

r L33 

X L3 + 
k 

<P 

l34 

X L4. 


B- 

-7 B5 

B- 

-8 B6 


BO B- 

•2 



+ 


L37 




L39 



B-4 

B2 

B2 

B-4 

BO 

B-2 

+ Y 

L31 

u 

CL1 , 
k 

+ Y 

L32 

u + 

CL2, 
k 

^31 

w 

L1 k 

BO 

B-2 

B-1 

B-1 

B-1 

B-1 

+ Y r 

LW 3 3 

w ^ 
L \ 

+ y T 

LW 34 

X 

+ Y t 

LW 35 

W „ 

L5, 


B-2 BO B-2 BO B-2 


+ Y 


Lw 


w 


36 


L6, 


+ ♦. 


L31 LI 


+ t 


L32 L2, 


+ 


L33 L3, 


B-2 

A. 

BO 

B-2 B-1 

B-1 BO 

B-2 B-1 

B-1 

X L4 

k+1 

<P 

L41 

x „ + <p 
LI, Y L42 
k 

X L2 + *1.43 
k 

X L3 * *L44 
k 

*L4. 


B-6 

B5 

B-7 

B6 

B1 

B-2 

+ <p 

L47 

X 

+ <j> 

L48 

x ^ + 
L8 v 

<P 

L49 

X L9, 


B-3 B2 B3 B-4 B-8 B-7 

+ Y L41 U CL1 + Y L42 U CL2 + Y L43 U CL3 
* k k 


B1 

B-2 

B1 

B-2 

BO 

B-1 

+ Y t 

LW 41 

w 

L1 k 

+ Y t 

LW 43 

w 

L3 k 

+ Y r 

LW 4 4 

w 

L4, 

i 

BO 

B-1 

B1 

B-2 

B1 

B-2 

+ \ 
LW 45 

w _ 
L5 k 

+ Y t 

Lv* 

46 

w „ 
L6 k 

+ ip 

L41 

x 
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B1 


B-2 


+ \h n 

L42 L2, 


B1 

+ 


B-2 


T\ 

L43 L3, 


BIO (DP) 
* 

B-1 

B-2 B-2 

B-1 B-1 

B-2 

B-2 

B-1 


d> 

LSI 

X L1 + ^L52 
k 

X L2 + *L53 
k 

x 

+ <J> 
y L54 

^4 


B1 0 (DP ) B-8 

B5 

B-9 

B6 

B-1 B-2 

+ X _ 

L5 k 

* *1.57 X L7 * 
k 

d> 

l58 

X* 

*1.59 X L9„ 
k 

B-5 

B2 

B-1 

B-2 

B-1 

B-2 

+ Y L51 

U CL1 

+ 

Lw r 

W L1 

+ Y Lw 

W L3 


k 

51 

k 

53 

k 


w 

to 

B-1 

B-2 

B-1 

B-1 

B-2 

Y r 

W M 

+ Y 

w _ 

+ Y 

w 

LW 54 

L4 k 

LW 55 

L5 k 

LW 56 

L6. 


B-1 B-2 B-1 B-2 B-1 B-2 


+ ♦ 


LSI 



+ t 


L52 



L53 



B-2 


+ 


L55 



B6 (DP) 

A 

B-1 

B-2 B-2 

B-1 B-1 

B-2 B-2 

B-1 

X L6 

k+1 

= <p 

L61 

X L1 + *L62 
k 

X L2 + ^L63 
k 

X L3 * *164 
k 

*L4. 


B6 ( DP ) 

B-8 

B5 

B-9 

B6 

B-1 

+ x _ 
L6 k 

+ & 
Y L67 

X L7 + ^L68 
k 

X + 

L8 k 

Y L69 

B-5 

B2 

B-1 

B-2 

B-1 

B-2 

+ Y L61 

U CL1 

k 

+ Y t w 

Lw Ll 

61 k 

+ Y t 

LW 63 

W L3 


65 


B-2 


B-2 
+ Y Lw 


64 


B-1 
W L4 


B-2 
+ T 


B-1 

w_ 


LW 65 L5 k 


B-1 
+ T LW 


W, 


66 


L6, 


B-1 

B-2 

B-1 B-2 

B-1 

B-2 

y L61 


+. ib n 

V L62 L2 k 

+ ^ 
V L63 

n L3. 


B-3 BO 


^ ^ n ^ 

L66 L6 


B5 

BO 

B5 

B-1 B6 B-2 B7 

A 

X L7 

k+1 

= <t> 
Y L77 


*1,78 X L8 k + Y L73 U CL3 k 

B-2 

B-2 

B2 

B-4 

^\+1 


+ Y L92 

U CL2 

k 


2. Winds 


B-2 

BO 

B-2 BO 

B-2 


= <p 

Lw n 

WL1 k + * LW 12 

W L2 


B-2 

BO 

B-2 

BO 

B-2 


= 4> 

LW 21 

W T 1 + 

<(> 

LW 22 

X 

B-2 

BO 

B-2 

BO 

B-2 

A 



<t> 

LW 32 


w „ 
L3 k+1 

= 

LW 31 

X* 

x 

B-1 

BO 

B-1 



w 

L4, , 

k+1 

*©■ 

11 

X 




B-2 


W, 


L3 


k 


66 


I 



B-1 

B-1 

a 


" L5 k + i 


B-2 

A 

B-2 

W L6 

k+1 

ll 

* 


B- 1 B1 B-2 


B-2 

B-2 

A 


w „ 

= w 

L7, 

L7 

k+1 


3* Biases 


B-2 (DP) 

A 

B- 2 ( DP ) 

A 

BO 

B-2 

B-6 

B4 B-6 B4 



s, 

V 

L \ 


V + f V 

L3. Lb L4 

k 14 k 


B-6 

B4 

BO 

B-2 

BO B-2 




B-1 B-1 BO B-2 

+ f . V + f V 

lb 18 L8 k Lb ,9 L \ 


B4 ( DP ) B4 ( DP ) B4 B-2 B-2 B4 B-2 B4 


L5 


k+1 


= b + f v + f v +f v 
L5 k Lb.., L1 k Lb 53 L3 k Lb 54 L4 k 


B-2 B4 B4 B-2 B4 B-2 

+ f v +f v +f v 
Lb cc L5 v Lb cc L6 Xb L7 

55 k 56 k 57 k 


B3 B-1 B4 B-2 

+ f , V + f V 

Lb 58 L8 k Lb 59 L9 k 
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L6 


= b , + f 


k+1 


61 


B-2 

B-6 

B4 B-6 

V 

+ f t 

V + f 

L \ 

Lb 

63 

L3 Lb 

k i 


64 “’k 


B-6 B4 BO B-2 BO B-2 

+ f . V 


TC + f V + f V 

Lb 65 L5 k lb 66 L \ Lb 67 L7 k 


B-1 B-1 BO B-2 

+ f , V + f V 

Lb 68 L8 k bb 69 L9 k 


B-2 (DP) B^2(DP) B-4 B-2 B-10 B4 B-10 B4 B-10 B4 

73 k 74 “ *k 


bb vr H * f 


bb 71 H + H 3 H + X, X * X 5 X 


B-4 B-2 B-4 B-2 B-5 B-1 

+ f , v + f V +f u 

Lb 76 L6 k “> 77 L7 k bb 78 L8 k 


B-4 B-2 

+ f V 

Lb 79 L \ 


B- 2 ( DP ) B-2 (DP) 

+■ ys 

B-4 

B- 

-2 

B-10 

B4 

B-10 

B4 

X., * H + 

fLb 91 

V + 

L1 k 

f Lb 

93 

v + 

L3 k 

f Lb 

94 

V 

L4 i 

B-10 

B4 


B-4 

B-2 

B-4 

B-2 


* X 5 

V 

L5 k 

+ 

fLb 96 

V 

L6 k 

+ 

97 

V 

L \ 


B-5 

B-1 


B-4 

B-2 




+ f T . 

V 

+ 

f 

V 




Lb 98 

L \ 


Lb^ 

99 

L9 k 
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OQ < 


4. Measurement Vector 


B-2(DP) BO B-2 B-1 B-1 


BO B-2 


B-2 (DP) 


L6 


k+1 


' °L*,. *L1. ... *• °L* X L2 * C Lx fi , X L3. , + V , 

62 k+1 63 k+1 k+1 


61 k+1 


B-2 (DP) BO B-2 

A ^ 

y. 


B-1 B-1 


BO B-2 


L V . ' ° LX 71 *“», + °**, 2 + Clx 73 H *, 


B-1 B-1 B-7 B5 


BO B-2 


+ + ^,7 X, X, 


BO B-2 


BO B-2 


B-1 B-1 


+ C W +C W +C w 

LW 71 L1 k+1 LW 73 L \+1 LW 74 L4 k+1 


B-1 B-1 


BO B-2 


B- 2 ( DP ) 


+ C w + c w _ + b _ 

Lw L5 Lw L6 , L7 . 

75 k+1 76 k+1 k+1 


B-1 (DP) B-1 (DP) B1 B-2 


^L8 X L2 + °Lw W L1 

k+1 k+1 81 k+1 


BO 
+ c 


B-1 

A 

W_ 


Lw ’ L4, „ 

84 k+1 


BO 

B-1 

A 

B1 

B-2 

A 

+ c 

w _ 

+ C 

w 

LW 85 

L5 k+1 

LW 86 

L6. 


k+1 


-2 (DP) BO B-2 


B-1 B-1 


BO B-2 


H *, " X , + X 2 X ^ k+ , + H *, 


BO 
+ c 


B-1 

A 

w 


Lw 85 L5 k + , 


B1 
+ c 


B-2 


L “86 X + 1 
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B- 2 ( DP ) BO B-2 

Ht, ' '“9, H*, 


B-1 

B-1 

A 

c 

Lx^ 

94 

X L4 

BO 

B-2 

A 

c 

LW 9! 

W 

LI , 

BO 

B-2 

CLW 96 



B-1 B-1 

BO 

B-2 

A 


A 

CL *92 *“*♦! 

+ C T 

LX 93 

X L3, 

B-7 B5 

A 

BO 

1 


+ C + X + c v 

LX 97 L7 k+ 1 *-99 H+1 


B- 1 B- 1 


B- 1 B- 1 


+ c w . + c w 

Lw L4 Lw„ r L5 

94 k+1 95 k+1 


B-2 (DP) 


+ b. 


L9, 


k+1 


B. Lateral Predictions 


1 . State 


BO 

A 

BO 

BO 

B1 

B-1 

B1 

B-1 

B-1 

B1 

V, 

= <t> 

Y 1 1 

\ + 

<f> 

Y 14 

v 

*15 x 5„ - 
k 

Y 1 1 

\ 


B1 

B-1 

B- 

2 B2 

BO 

BO 




- C 

13 ^3 

k 

- ^ c 

14 4 

k 

- \p £ 

15 5, 

k 



B2 

B2 

BO 

B2 

B3 

B-1 

B3 

B-1 



= 

r 2 1 


\ 

+ <p 
y 24 

X 4 + *25 
k 

\ 



B1 

B1 

B3 

B-1 

B3 

B-1 

BO 

B2 


\ ■ *» \ ■ *» \ - * 24 v» 
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BO 

BO BO 

BO 

BO 

B1 

B-1 

B1 B-1 B1 B-1 

X 3 

k+1 

= *31 X 1 

+ <b 

33 

k 

\ 

y 34 


*35 X 5 + *37 X 7 
k k 


B-2 

B2 

B-2 

B2 

B1 

B-1 B2 B-2 


+ Y 31 


Y 32 

U C2 

^k 

+ Y 

W 31 

w + Y w 

'k W 33 3 k 


B2 

B-2 

B1 

B-1 

B-1 

B1 B1 B-1 


+ Y W w 4 + Y w w 5 " *31 S " *33 S 
34 k 35 k k k 


B-2 B2 BO BO 


I* C. - 

34 4, 35 5, 

k k 


B-1 

* 

x. 


k+1 


B-1 BO 

B-1 

BO 

BO 

B-1 

BO 

B— 1 

BO 

B-1 

*41 X 1 

+ d> 

k 43 

\ 

+ *44 


+ *45 

\ 

+ <)> 
*47 

\ 

B-1 

B2 

B-3 

B2 

BO 

B-1 

B1 

B-2 


+ Y 41 


Y 42 

U C2 + 
ux k 

Y 

w 

41 

\ 

+ Y 

w 

W 3 

43 k 


B1 

B-2 

BO 

B-1 

B-2 

B1 

BO 

B-1 


+ Y 

W 44 

\ + 

Y 

W 45 


ip 

r 41 

v 

*43 




B-3 B2 B-1 BO 

C, " ¥. c C c 

44 4, 45 5, 

k k 


B-1 B-1 BO B-1 BO BO B-1 BO B-1 BO B-1 

X = * x + 4> x + <t> X + * X + ♦ x 

k+l k k k k k 


B-3 B2 B-3 B2 BO B-1 BO B-1 


+ Y 51 U C1 + Y 52 U C2 + V V + Y w„ W 3 
k k 51 k 53 k 


71 



B1 B-2 BO B-1 B-2 B1 BO B-1 

+ Y w + y w C -ip C 

“54 \ ”55 \ 51 V 53 3 k 


B- 3 B2 B-1 BO 




k+l 


BO 

BO B- 2 B2 

BO 

BO 

B1 B-1 

B1 B-1 

<t> 

61 

V + *2„ 

k k 

+ <f> 

63 

\ 

+ <b x 

*64 4, 

lc 

+ d> x 

65 5, 

lr 


B6 ( DP ) B1 B-1 B-2 B2 B-2 B2 

+ X + y u , + Y u 

6 k 67 7 k 61 cl k 62 c2 } 


B1 B-1 B2 B-2 B2 B-2 B1 B-1 

Y w + Y w +Y w +Y w 

W 61 \ W 63 3 k W 64 4 k W 65 5 k 


B-1 B1 B1 B-1 B1 B-1 B-2 B2 

*61 s - * 62 s - * 63 s - * 64 ,t 4 . 


BO BO B-2 B2 

^65 S ~ ^66 C 6 
k k 


B-1 
x. 


B-1 B-3 B2 


7 = x + Y.,_ u „ 

? k+1 7 k 72 C2 k 


2. Winds 


B- 

A 

•1 

BO 

B-1 

B1 

B-2 

W 


= <t> 

W„ 

+ <t> 

w 

1 

k+l 

W 11 

\ 

W 12 

2 k 


B-2 

B-1 

B-1 

BO 

B-2 

w o 

= <t> 


+ 4> 

w 

2 k+i 

w 

21 

V 

w 

22 

2 k 
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B-2 

B-1 

B-1 

BO B-2 BO 

B-2 

A. 





W 

k+1 

w 

31 

\ 

+ <fr + <f 

w 2 w 

32 k 33 

\ 

B-2 

BO 

B-2 



w 

k+1 

4> 

W 44 

\ 



B-1 

B-1 

B1 

B-2 


* 





W 5 

k+1 

\ 

+ <|> w 

w 6 

56 k 


B-2 

BO 

B-2 



w 

k+1 

* 

W 66 

\ 



Biases 





BO ( DP ) 
* 

BO (DP) BO 

A 

BO 


b * = 
y k+1 

b. 

+ g * 

V 


B2 ( DP ) 

B2 (DP) BO 

B2 


A 

A 




b i 

*k+1 


+ g * 

V 

2 k 


B-2 ( DP ) 

B-2 (DP) 

B1 B-3 


y k+1 

= b f 

+ 

9f \ 
y k 



4. Measurement Vector 


B-3(DP) B-2 BO B-1 B-1 B-1 B-1 

A A A A 

V = c + c x^, + c x_ 

6 k+1 x 63 3 k+1 x 64 4 k+1 X 65 5 k+1 


B-1 B-1 B-1 B-1 BO B-2 

A A A 

+ C + C V/ + c w 

x 67 ? k+1 W 61 1 k+1 W fi3 \+1 
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BO B-2 B-1 B- 1 B-2( DP) 

ys A 

+ c w_ + c w_ + b 

W 64 4 k+1 W 65 5 k+1 fy k+1 


C« Euler Prediction Matrices 


sin 


0 = 


sin 

cos 


B2 B-2 

B-1 

A 


9 + x 

1 

o L.1 


k+1 

TT 


BO B1 


sin 

cos 


sm . 
cos \ X 1 


k+1 


sin 

cos 


♦ = 


B2 B1 

f 

sm 

cos V 2. „ tt 

k+1 


ES22 = sin <f> sin 9 sin ^ + cos <j> cos ^ 


EB31 = -sin 0 


EB32 = sin 4> cos 9 


1 EB3 3 = COS ^ COS U 


1 ~ 1 C + 1 q 

ES3 1 EB31 a EB33 a 

o o 


1 =1 ( 1 i p 

ES3 3 EB31 v a ; EB3 3 a 

o o 


74 


l 



D • Altitude and Path Parameter Prediction 
1 • Altitude Prediction 


If not GSENG, Go to 2 . 


BO 

A 

n. 


BO 


BO 


B— 1 


= (1 ^ + Sot ) x ^ 

L6. 4 ES31 o L2 „ 

k+1 k+1 


BO 

A 

+ 1 


ES3 1 


B-2 

BO BO 

BO 

B-2 

BO BO 

A 

A 


A 

A A 

X 

CY + (1 - 

CY ) 

x ^ 

+ 1 x_ 

LI. „ 
k+1 

o ES33 

o 

L3 n . 
k+1 

ES32 3, , 

k+1 


B13 B6 ( DP) B8 


H = 



U 

o 


2* Sideslip Parameters 


If not LOCE, skip over 3 and £ computations 

K K 


B-1 

A 

6 . 


B-1 

A 


k+1 


k+1 


aBO 
1 


ES2 2 


B-1 B-1 


k+1 k+1 


B-1 


= 3, - w_ 


B-1 

A 

- 


k+1 


k+1 


7 5 



V. INITIALIZATION 


A. Computation of a , 9 
o o 

u ok ~ calibrated airspeed selected by pilot (knots) 
WEITO = empty fuel weight (LBS/10) 

B7 

WEIGHT 

7000 ~" = WEIT0 / 100 + TPQ/1 000 
If (U ok 4 127) ' then 


B-6 

B-14 

B8 

B-6 (DP) 

a i = 

0.000056791 

(U J 

ok 

- 0.007992711 

B1 

B-7 

B8 

B1 (DP) 

a 2 = 

-0. 003278089 

(U 

ok 

+ 0.465299649 

B1 

B-7 

B8 

B1 (DP) 

b 1 = 

-0.003272937 

(U 

ok 

+ 0.559452617 

B8 

B0 

B8 

B8 ( DP ) 


b = 0.061661686 (U , ) - 18.39413409 
^ ok 

Otherwise 

B-6 B-14 B8 B- 6 ( DP ) 

a = .000029097 (U , ) + -0.004475519 
1 ok 


B1 

B-7 

B8 

B1 (DP) 

a 2 = 

-0.001 527618 

(U J 

ok 

+ 0.242989939 

B1 

B-7 

B8 

BI (DP) 

b i = 

-0.004484078 

(U J 

ok 

+ 0.713267456 


76 



B8 BO B8 B8 ( DP ) 

h> 2 = 0.195679600 (U Qk ) - 35.41440925 


B1 


B-6 B7 


B1 


.WEIGHT, 

a 3 ■ a i ‘-iooo-’ + a : 


B8 B1 


B7 


. .WEIGHT. 
b 3 - b , ( -7ooo-> * b 


B8 

2 


Y = GAMO where GAMO is entered value 
o 


a 

o 


(a Y + b ) 
3 o 3 


0 = ct + y 

o o o 


BO 

S = sin ot 
a 

o 

BO 

C = cos ot 
a o 

o 


BO 

SY = sin Y 
o o 


BO 

CY = cos Y 
o o 


BO BO 

c (6,1) = -cos Y 
Lx o 


B-1 BO 

c (6,2) = -sin Y 
Lx o 


77 



BO BO 

c (6,3) = cos Y 

LX O 


B1 

SCY = 1/CY 
o o 


B1 

T - tan Y 
Y o 

o 

BO 

S n = sin 0 
y o 

o 


BO 

C n - cos 0 
U o 

o 

B8 B7 

U o = U o (1 - 687 8) 


B-7 B1 B8 
UOI = 1 ,/u 


B, Computation of T , Th (Routine 

o o 


B-7 


B-1 4 


B7 


a = 0.000034840605 1 T ,) 

1 1000 ' 


B1 


B-6 


B7 


a = .01299031000 (-— + 
2 1000 


B-1 B-8 B7 

b = -.00245252350 (WglgHT) + 
1 1000 ' 


7B 


INIAT) 

B-7 

.00201123645 

B1 

29543500 

B-1 

.2791 23915 


I 



B7 


BO 


B7 


B7 


fc> 2 = .405932250 (^|~^) - 29.95779350 


B1 B-7 B8 B1 


a., = a, U , + a„ 
3 1 ok 2 


B7 B-1 B8 B7 

b, = h U , + b„ 

3 1 ok 2 


B5 B1 B6 B7 


T = a Y + b, 
o 3 o 3 


B6 B1 


1000 lbs 

where Y is in deg 
o 


B6 


Th = 1.51994 T + 5.257256565 
o o 


C. Glides lope Parameters 


Z INT “ 0 


J GSC “ 0 


k fl = 0 


e LX6 k = 0-125 


V 0 


b A = 0 
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0. Flare Parameters 


EZ5 = 

0 


EZ7 = 

0 


J Flare 

= 0 


> 

< 

* 

II 

0 


k fl “ 

0 


1! 

O 

>- 

r 

y 

0 


B- 1 

B6 

B-7 

II 

$ 

-25 

UOI 

B-2 

B5 

B-7 

Av 

FR 

1 .25 

UOI 

B- 1 ( DP ) 

B-3 

B-7 

Av 

FD 

( 1 .25At) UOI 

B6 

B1 3 

B-7 

h 

FLR 

HFLR 

F 

UOI 

B6 

B1 3 

B-7 

x = 

FLR 

XFLR_ 

F 

UOI 


B6 

BO B6 

B6 

h i 

h 

FLR 

- T X 

y FLR 

0 

V2 p 

BO 

BO 


T - T 

yTD y 

0 
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B6 


B6 


Ax f = h i 


/2 P 


B6 

x 


FLR 


B6 ( DP ) BO B6 


h _ = T 
f Y_ 


(x 


FLR 


B6 
- 2h 




BIO(DP) B6 

X = 

FLR FLR 


h f • 


Z L6 


h 

h 

h 


Lx3 

Lx4 

Lx6 


= h 

= 0 

(i) 

(i) 

(i) 


h LxOi3 

h Lx0i4 

= h T 
Lxi6 


> 


for i = 1, 2, 3 


E. Wind Parameters 

B-1 B8 B-9 

DT( 4 ) = U * .00099758 

o 


l> . = -sin Y <1' . + cos Y 

L5i o LI i o L2i 


+sin Y . 

o L3i 


|> _ . = -cos Y <l> . - sin Y . 

L6i o LI l o L2i 


+cos Y . 

o L3i 


for i = I, 4 and i =7,9 
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VI 


PATH CAPTURE CRITERIA AND LOGIC 


A. Pre-Capture Eq uations 


B-2 BO 

® = Asin (sin 9) 
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APPENDIX B 


MLS SIGNAL PROCESSING 

The Microwave Landing System (MLS) used in the flight tests provided three 
signals from which the aircraft's position, relative to the MLS installation, could 
be determined. Two signals were from the MLS angle receiver-azimuth angle relative 
to runway centerline (+60 degrees coverage) updated at 13,3 Hz and elevation angle 
relative to the local level at the elevation antenna site (1 to 20 degrees coverage) 
updated at 40 Hz. The third signal was from the MLS DME receiver representing the 
distance or range to the DME antenna site and was updated at approximately 20 Hz. 

The location of these MLS ground antenna sites is shown in figure B1 . Using these 
signals, the aircraft position in x-, y-, and z-coordinates relative to an origin at 
the MLS azimuth antenna phase center was computed as follows (see fig. B1 ) : 

The range, R, to the origin was computed using the equation 


R 


R 

m 


y_. ir( sin AZ + R 
DME o 


where R was the MLS DME receiver measurement, y the DME ground antenna 

y-coordinate, AZ the MLS azimuth measurement, and R a constant accounting for the 
x-coordinate of the DME ground antenna, x D , anc * cable length bias, R c between 
the aircraft's antenna and the MLS DME receiver. R is given as 


R 

o 


X DME 


R 

c 


The equations for the x- , y-, and z-coordinates are 


x = g (g 2 - h)^ 2 
y = -R sin AZ 


Z E = Z e£ + tan E£ / (X " X E* )2 + (y ‘ y E* )2 


where x and Z E £ are coor< ^ nates the MLS 9 round elevation antenna's 

phase center aria g and h are defined as 


2 

9 = * EJl sin El 

2 2 2 2 2 2 2 
h = (x ^ + sin + y ~ R cos E& - 2y y^ sin e£ 


where E £ is the elevation angle measurement. 
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When the aircraft approaches the touchdown point on the runway, the elevation 

signal from the MLS becomes invalid due to the aircraft flying out of the elevation 

coverage (the antenna siting is usually 250 to 400 feet from the runway centerline 

and BOO to 1000 feet from the threshold; a nominal touchdown point is on the order of 

1500 feet). Thus, the altitude information provided by the elevation signal was 

replaced by an alternate measurement, radar altitude (h ). The x-, y-, and 

z-coordinates were computed when using h as R 

R 

Pitch Runway Slope 

Correction Correction 

a 

z = h + h - 0.380 - 0.0023 (X - x ) 
n R R) o 


X 


'R (1 


. 2 
sm 


AZ) 


2 

z 


y - -R sin AZ 


where h R ^ is a constant to transform the radar ^altitude measurement to the roof 
antenna, 0 is the pitch attitude in degrees, X is an estimate of the x-coordinate 
from the third-order complementary filter to be described later and x is a con- 
stant. The values used for h ^ and x were, respectively, 8.86 and 8049.5 feet. 
If x was less than the runway correction was not used. 

The switch from El to h R occurred when the aircraft crossed over the runway 
threshold according to 

✓s 

if X < 9630 feet (approx. 400 feet from threshold) use h . 

R 

When the above logic became true, the z-coordinate computed from h was entered 
into the calculations through an easy-on with a 2 second turn-on time. The equation 
for the z-coordinate transition is 

z = 01 z, + (1 - <*) z 

hR E 


where a goes from 0 to 1 during the transition, z is the z-coordinate computed 

from h R , and z is the z-coordinate computed from Eft. 

The above equations will give the x- , y- , and z-coordinates values of the air- 
craft’s receiving antenna in the MLS coordinate frame. In general, control of the 
aircraft's center-of-gravity or c.g. position relative to some desired flight is 
desired. Thus, given the position of the aircraft’s antenna the coordinates of the 
aircraft's c.g. position were determined as 


_ 


— — 


— — 

X 


X 


X 

eg 




a 

y 


y 

+ T _ 

y 

eg 


MB 


z 


z 


z 

L cc ? J 


— — 


L_ a _ J 


1 04 


i 



where i.s an Killer transformation matrix which transforms a vector in the air- 
craft's body coordinates to the MLS coord i nates and x a , y^, and z are the body 

axis coordinates of the aircraft's MLS receiving antenna. The aircraft's MLS antenna 
was mounted above the cockpit on the roof of the aircraft resulting in coordinates 
defined as x^ - -35.05, y a = 0.83, = 6.07 (the positive direction of the 

antenna vector was from the aircraft antenna to the c.g.). T is computed using 
the attitude measurements — pitch, roll, and true heading (yawf^ The coordinates, 

X cg' ^cg' anc * Z cg' were ^ n P ut directly to the DIALS innovation equations . 

Th e x C g' ^cg' anc * Z cg coor< 3inates were also input to three third-order 
complementary filters — all identical and one for each coordinate axis. The form of 
these filters is shown in figure B2. The filters were developed for MLS processing 
(ref. 9) and were a part of the baseline software of the test vehicle. The filters 
provided estimates of the aircraft's position agd velocity in the MLS coordinate 
frame. The position coordinate estimates are X, Y, and Z. The velocity estimate 

of the MLS z-coordinate axis was designated as ^LS w hi- c h replaced the barometric 
altitude rate in DIALS during the flight tests. 
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Figure Bl. - WFC MLS Antenna Locations on Runway 22. 




1 07 


Figure B2. - Third Order Complementary Filter. 





APPENDIX C 


ALPHA BETA FILTERS 


Software filters were used in the flight control computers to provide filtering 
for the three body-mounted accelerometer measurements, in addition to that provided 
by the standard hardware input filters. These filters, activated on one run of the 
ight test data presented, were executed each minor frame or data input cycle of the 
flight control computer, that is, every 6.144 msec. The same filter was used for 
each of the three accelerometer measurements and thus the equations for each filter 
are the same. The equations for this filter are given as 


x = (1 - a ) x + y x 
n l Pn n 


X = X + 3 / (x 

n n-1 f / T n 


Pn 


X Pn+1 


x + T x 
n n 


where x^ is the measurement, x is the estimate or 
estimate of the rate of change of, the measured value. 

setting x equal to x and x equal to zero, 

constants were n n 1 


filtered value, and x is an 
The filter is initialized by 

The values of the filter 


a f = 0.13 
3 f = 0.009037 
T = 0.006144 sec 
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Figure 1. - The DIALS Block Diagram, 





Figure 2. - Small Twin-jet Test Aircraft. 
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Figure 3. - DIALS Horizontal and Vertical Profiles 
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Figure 4. Continued 
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Figure 8. 3 Glideslope Capture, Track, and Flare 





< 5 * 

deg 



Time, sec 


(b) Path and Control Variables 


Figure 8. Concluded 


121 







60 


80 

Timi 


(b) Path and ( 


Figure 9. 






0 
















10 60 


80 100 120 1 

Time, sec 


slope Capture for Aircraft Cl 


1 46.2 r 


V 0 

knots 


126.2 

106.2 

8 



deg/sec 
9 

deg 


- 





deg 


fps' 



fps' 



0 20 40 60 80 100 120 140 160 180 200 

Time, sec 


(a) Measurement Variables 


Figure 14. 4.5° Glideslope Capture, Track, & Flare for Mild Winds 
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